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Abstract. We define a Cartesian product G x iJ for finite sim¬ 
ple graphs which satisfies the Kiinneth formula H^{G x H) = 
®i+j=kH'-{G) ® W{G) and so pgxh{x) = PG{x)PH{y) for the 
Poincare polynomial pg{x) = X]/c=o x{G x 

H) = x{G)x{H) for the Euler characteristic x(G) =Pg(~ 1)- The 
graph Gi = G X ATi is homotopic to G, has a digraph structure 
and satisfies the inequality dim(Gi) > dim(G) and G\. Hodge 
theory leads to the Kiinneth identity using the product fg of har¬ 
monic forms of G and H. A discrete de Rham cohomology and 
“partial derivatives” emerge on the product graphs. We show that 
de Rham cohomology is equivalent to graph cohomology by con¬ 
structing a chain homotopy. The dimension relation dim(G x 
H) = dim(G) -I- dim(iL) holds point-wise dim(G x H){x,y) = 
dim(Gi)(a:)-|-dim(iLi)(y) and implies the inequality dim(Gx H)> 
dim(G) + dim(iJ), mirroring a Hausdorff dimension inequality di¬ 
mension in the continuum. The chromatic number c(Gi) of Gi 
is smaller or equal than c(G) and c(G x H) < c{G) + c{H) — 1. 
Indeed, c(G x H) is the maximal n for which there is a A„ sub¬ 
graph oi G X H. The automorphism group oi G x H contains 
Aut(G) X Aut(iL). If G ~ iL and U ^ V are homotopic, then 
G X U and H x V are homotopic, leading to a product on ho¬ 
motopy classes. If G is fc-dimensional geometric meaning that all 
unit spheres S{x) in G are (fc — l)-discrete spheres, then Gi is k- 
dimensional geometric. And if H is /-dimensional geometric, then 
G X H is geometric of dimension (fc -|- /). Because the product 
writes a graph as a polynomial Jg of n variables for which the 
Euler polynomial e{x) = is eG{x) = /^(a;, ...,x) and 

x(G) = eG(—1), the product extends to a ring of chains which 
unlike graphs is closed under the boundary operation S defining 
the exterior derivative df{x) = f{5x) and closed under quotients 
G/A with A C Aut(G). By gluing graphs, joins or hbre bundles 
are defined with the same features as in the continuum, allowing 
to build isomorphism classes of bundles. 
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Figure 1. The product G x H oi the house graph G of 
dimension 1.4666 ... and a second graph H of dimension 
1.133... which is also homotopic to a circle produces a 
product graph of dimension 2.70238... which is homo¬ 
topic to a 2-torus. Like Hausdorff dimension in the con¬ 
tinuum, the dimension of the product is larger or equal 
than the sum 2.6 of the dimension of the two factors. The 
graph G X if is homotopic to a torus, has Poincare poly¬ 
nomial 1 + 2x + which is the product of the Poincare 
polynomials = 1 + x = Ph{x) of the factors. The 

Kiinneth theorem tells that if /, g are harmonic 1-forms 
representing a nontrivial cohomology class in H^{G) or 
fi^(if) respectively, then f{x) * 1,1 * g{y) can be used 
to construct a basis for H^{G x H) and f{x)g{y) can be 
used to build a 2-form spanning the 1-dimensional space 
H^{G X H). The discrete de Rham theorem tells how to 
get from the de Rham picture to the simplicial cohomol¬ 
ogy picture of G x if. While the linear space of 1-forms 
on G and if is 6 or 21-dimensional respectively, the de 


Rham 1-forms in f2^(G) (g)fi°(if) ©fi°(G) (8)f2^(if) build 
a6* 20-1-5* 21 = 225 dimensional space, while f2^(Gxif) 
has dimension 2196, which is almost 10 times more. And 

fi2(G) © n%H) © n^{G) © n\H) © 12^^) © n\H) is 
6*1 + 5*1 + 6*21 = 137-dimensional while 12^(G x if) 
has dimension 2880, which is the number of triangles 
in G X if and more than 20 times the dimension 137 
in the de Rham case. The discrete Eilenberg-Zilber 
theorem (10) assures that if^(G x if) is isomorphic to 
H\G)(^IF(h)®H%G)(^H\H) and that H\GxH) is 
isomorphic to if^(G) ©if^(if) here as if^(if) and H‘^{G) 
are 0-dimensional. The chromatic numbers of G and if 
are 3. The chromatic number of G x if is 5 as it contains 
a 4-dimensional clique if 5 and Theorem (24). 
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1. Introduction 

As acknowledged by nomenclature, Descartes concept of coordinates 
depends on the notion of a Cartesian product. Omnipresent in mathe¬ 
matics, it is already used in basic arithmetic to build number systems 
like the held of complex numbers or to dehne exterior bundles on man¬ 
ifolds. The Cartesian product allows to build and access higher dimen¬ 
sional features of geometric spaces. Many constructions in topology 
like suspensions, joins, hbre bundles, de Rham cohomology or homo- 
topy deformations would not work without the concept of a Cartesian 
product. Of course, we would like to have a product in graph theory 
which shares the properties from the continuum. The new graph prod¬ 
uct will achieve that. It will allow to use “coordinates” similarly as 
they are used in the continuum. For two arbitrary networks G,H - 
that is G and H are hnite simple graphs - the coordinates of the prod¬ 
uct graph GxH consists of all pairs of complete subgraphs of G and H. 
The exterior derivatives in G and H will play the role of the “partial 
derivatives” in the product and allow to build an exterior de Rham de¬ 
rivative on the product graph GxH. The actual exterior derivative on 
the product graph GxH operates on a much larger complex. The later 
is called the Whitney complex and is dehned by the simplices in the 
product. Having a product allows to work with “rectangular boxes”s 
in the product space rather than with simplices. Figure Q illustrates 
this for a small example. To take a picture from the continuum: the 
curl in the plane is an inhnitesimal line integral along a rectangle and 
uses the Leibnitz rule to relate the exterior derivative of the factors 
with the exterior derivative dfg — fdg of the product. A simplicial 
point of view of cohomology integrates around inhnitesimal triangles 
to get the curl, ignoring the product structure. De Rham establishes 
equivalence of the two pictures on any smooth manifold. In order to 
prove the Kiinneth formula, we will need to emulate the de Rham the¬ 
orem combinatorially. To do so, we explicitly construct a fc-form on 
the Whitney complex of Gx H from a fc-form on the de Rham complex 
dehned by the two graphs G, H. The chain map is concrete as we can 
from this construct explicit cohomology classes of the product from the 
cohomology classes in G or H. The relation is what one calls a chain 
homotopy. The de Rham picture will be useful if we work with “dis¬ 
crete n-manifolds” obtained by gluing together local charts of products 
Ui = Gil X ■ ■ ■ X Gin of networks or when working with “hbre bundles” 
IT ■. E ^ M obtained by gluing locally trivial charts UiX G of networks 
above a discrete manifold covered with charts Ui. The automorphism 
group of the hbre G will then play the role of a gauge group on E, as 
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it does in the continuum. 

When looking at graphs as geometric structures without dimension re¬ 
striction, an amazing similarity with the continuum emerges. It turns 
out that in the discrete, one can access the local structure of space G 
directly as the complete subgraphs of G. These simplices can serve as 
fundamental entities playing the role of “points”. Such insight has been 
promoted already in [7] but most of the time still graphs are treated 
as one dimensional simplicial complexes. An illustration on how the 
change of view point allows to emulate results from the continuum is 
the hxed point theorem of Brouwer and Lefschetz which looks identical 
to the result in the continuum I2D1. Many concepts become elementary: 
cohomology is part of hnite dimensional linear algebra, to compute val¬ 
uations, generalized volumes, one needs integral geometric tools in the 
continuum, while in the discrete is is just count of complete subgraphs. 
The discrete Hadwiger theorem [12] is much easier than the continuum 
version: the numbers Vi{G) of ^-dimensional simplices is a basis for the 
linear space of valuations. Differential forms are just functions on a 
simplex graph and Stokes theorem in the Whitney complex is a tau¬ 
tology; it becomes only less obvious when Stokes is considered in a de 
Rham setup. Much Intuition about higher dimensions can be obtained 
inductively. There are notions of dimension, cohomology, homotopy, 
cobordisms, ramihed covers, degree and index, spheres, geodesic lines 
and curvature which lead to results mirroring the results in the contin¬ 
uum. This is not only true for nice geometric graphs but for general 
undirected networks or hnite simple graphs - and much works with¬ 
out any exceptions. The notion of homotopy of graphs for example 
immediately leads to the homotopy of a graph embedded in an other 
graph and so to homotopy groups for general hnite simple graphs. The 
product dehned here will actually help to dehne the homotopy groups 
as graphs might be too small at hrst to have spheres embedded, so 
that the graph should hrst be rehned. But the Hurewicz homomor- 
phisms from the homotopy groups 7ik{G) to the cohomology groups 
H^{G) are then so explicit that one can even watch it happen: just 
apply the heat how to a /c-forms with support on the fc-simplices on 
the embedded fc-sphere. It converges to a harmonic form which by 
Hodge theory represents a cohomology classes. In the continuum, such 
a proof requires de Rham currents, generalized diherential forms which 
require some functional analysis. In the graph case, the heat how is 
just a linear ordinary diherential equation of the type studied in in¬ 
troductory linear algebra courses. The dehnition of homotopy groups 
which are relevant in coloring questions for graphs ^21 |2H] allows to 
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work with spheres in graph theory in the same way as in the continuum. 

The Euclidean product is not only essential for dehning fundamental 
objects like hbre bundles, it is also needed to construct spaces which 
have the same properties than classical manifolds. Examples of such 
properties are dimension, homotopy, cohomology or Euler characteris¬ 
tic. The goal of this note is to give such a product, allowing the use 
tools like discrete hbre bundles in graph theory. We will see that this 
can be done purely algebraically: as we can glue graphs together, this 
gluing carries over to the product allowing to build discrete bundles. 
And if the hbres carry an automorphism group A, we get discrete ana¬ 
logues of principle bundles on which an enlarged gauge group acts. 

A Cartesian product for a category of a geometry needs to be dimension- 
additive, it needs to induce a product on the homotopy classes, it needs 
to be Euler characteristic multiplicative, it must satisfy the Kiinneth 
formula equating the tensor product of the cohomology rings with the 
cohomology ring of the product and it must have the property that the 
automorphism group contains the automorphism groups of the factors. 

As for graphs, no previously dehned product shares these properties. 
The standard Cartesian product “ x ” of the cyclic graph C 4 with C 4 
for example is a graph of dimension 1. Its vertices are the Cartesian 
product of the vertices and two points {x,y), {u,v) are connected, if 
{x, u) E E OT {y, v) G E. The cohomology of the standard Carte¬ 
sian product has little to do with the cohomology of the factors: the 
Betti numbers of 6 * 4 “ x ’'C 4 for example is {bQ,bi) = (1,17) while the 
Betti number of our product is (1, 2,1), which is identical to the one of 
the two-dimensional torus in classical topology. The dimension of the 
traditional product is 1 while the dimension of our product is 2 . 

Also other constructions like the “tensor product”, the “direct product” 
or the “strong product” of graphs do not have the topological properties 
we want. The reason for the shortcomings of all these products is that 
they do not tap into the lower dimensional building blocks of space: 
these are the simplices = complete subgraphs in the graph case. What 
in the continuum has to be done with sheaf theoretical constructs, is 
already pre-wired in the graph as we can access the lower dimensional 
simplices as “points”. The new product has the property for example 
that C4 X Cs X Cii is a 3-dimensional torus of Euler characteristic 0 
and Betti vector (1, 3, 3,1). It is a triangularization of and each unit 
sphere S{x) is a 2-dimensional sphere of Euler characteristic 2. The 
product has practical use as it allows to construct high dimensional 
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Figure 2. The Mobius bundle is an example of a non¬ 
trivial bundle. Locally, it is a Cartesian product of Pi x 
where Pk is the 1-dimensional line graph with k + 1 ele¬ 
ments. But the graph itself is not orientable. Its coho¬ 
mology has the Betti vector b = (1,1, 0), like the circle. 
While not distinguishable from the cylinder by cohomol¬ 
ogy nor homotopy nor dimension, its topology is different 
as its boundary (the subgraph generated by the vertices 
for which the unit sphere is not a sphere) is not connected 
and it is not orientable. As in the continuum, a bundle 
can be constructed by taking a base graph G, cover it 
with a nice open cover such that its nerve is homotopic 
to G. Now build the product graphs and make sure that 
the transition maps are homeomorphisms in the sense of 

pj]. 


geometric spaces from smaller dimensional ones. The product works for 
any pair of graphs and leaves geometric d-dimensional graphs invariant, 
graphs for which all unit spheres are (d — l)-dimensional homotopy 
spheres dehned in [2B]- The product will again have this property. 
Cohomology, dimension and homotopy properties of the product are 
identical to the properties in the continuum. Even in the case of fractal 
dimension, the dimension formula matches the corresponding product 
formula in the continuum [6] (formula 7.2) for the Hausdorff dimension 
of arbitrary sets in Euclidean space. There are more analogues: for 
Hausdorff dimension, there are sets of dimension zero for which the 
product has dimension 1. While graphs of dimension 0 are geometric, 
so that one gets equality, there are sequences of graphs like = 
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Figure 3. The product K 2 x K 2 is the wheel graph 
Wg. Its boundary is a 1-sphere. The graph K2 x K2X K2 
is a 3-dimensional ball with 27 vertices and 98 edges. It 
resembles a cube as it is 3-dimensional and has six 2- 
dimensional faces of the form K 2 x K 2 . The product 
{K 2 X K 2 ) X ii '2 is a 3-dimensional ball with 99 vertices 
and 466 edges. Its boundary is a 2-sphere. As {K 2 x 
K 2 ) has been obtained by returning back to geometry 
before computing the new product {K2 x K2) x K2 is a. 
rehnement of K2 x K2 x K2. similarly as G x is a 
rehnement of G. The boundary of K2 x K2 x K2 by the 
way agrees with the enhancement O x Ki oi the octahe¬ 
dron. 

XqXi + Xq + Xi + - ■ --^Xn for which dim(G„) —)■ 0 and dim(G„ x G„) —)■ 1, 

mirroring the continuum again. 
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An other useful feature of the product is that it allows to rehne graphs: 
the enhanced graph G x A'l is a barycentric rehnement of G and this 
can be repeated: the sequence G x Ki,{G x Ki) x Ki, - ■ ■ produces 
a hner and hner mesh whose dimension converges to an integer and 
honor the original symmetries if A C Aut(G) is a subgroup of the au¬ 
tomorphism group then A still acts on the rehnement. The later is 
important when looking at discrete principal bundles. Moreover, while 
G/A is in general no more a graph, the quotient G x K^/A is, so that 
we can as in the continuum form covering spaces ramihed over rather 
general subgraphs. For example, while A = Zn acting on Gn has a 
quotient Gn/A which is only a chain and no more a hnite simple graph, 
the quotient (G„ x Ki x Ki)/A is a graph. As an other example, the 
quotient 0/A of the octahedron modulo the Z 2 action given by the an¬ 
tipodal rehection group A = Z 2 is no more a geometric graph without 
boundary as it is the wheel graph IF 4 for which we can see the octahe¬ 
dron as a double cover ramihed over the equator (Riemann-Hurwitz in 
the discrete is just the Burnside lemma considered simultaneously for 
the various simplex sceletons as noted in [30]), but {O x Ki x Ki)/A 
is a geometric graph, a discrete projective plane. As unit spheres of 
geometric graphs are spheres, the product can be used to construct 
new spheres. By taking rehnements and then taking quotients, one 
can get more general discrete graphs like the projective plane in the 
simplest case. The product construction can also help for other con¬ 
structions, like constructing joins or building discrete Hopf hbrations 
in arbitrary dimensions. As every unit sphere in the product has the 
structure 6 Bi x B 2 UB 1 X 6 B 2 which is the union of two solid tori glued 
along the torus 6 Bi x 6 B 2 . In the case when taking the product of 2- 
dimensional geometric graphs, we get so 3-dimensional spheres which 
have naturally the same Hopf hbration structure as in the continuum. 

Classically, the notion of a homotopy of two continuous maps /, g : 
X —)■ R is dehned using the product: if there is a continuous map 
F : X X [0,1] —)■ F, such that F{x, 0) = f{x) and F{x, 1) = g{x), then 
/, g are called homotopic. This could be done now also for graphs: two 
graph homomorphisms f,g:G^H are homotopic, if there exists a 
line graph and a graph homomorphism F from G x Ln to H such 
that F{x,0) = f{x) and F{x,n) = g{x). We have suggested in [25] 
an other dehnition: two graph homomorphisms f,g:G^H are ho¬ 
motopic if the “graph of the graph homomorphisms” are homotopic 
graphs. These graphs of homomorphisms have as vertices the union 
V{G) U V{FI) and as edges all pairs {xi,X 2 ) G E{G), {yi,y 2 ) G E{H) 
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and pairs {x, y) with y = f{x). Two graph homomorphisms are now ho- 
motopic, if the corresponding graphs are homotopic as graphs. We be¬ 
lieve that these two dehnitions are equivalent, but have not yet proven 
this. In any case, we have the Whiteheads theorem that if there is 
a graph homomorphism / : G —)■ if which induces isomorphisms on 
homotopy groups 7r„(G) —)■ 7r„(ii), then G,H are homotopic. And as 
in the continuum, having isomorphic homotopy groups does not force 
a homotopy equivalence. Now, with a product we can take the same 
standard counter example G = x and H = x S^, where P^ 
are graph implementations of the i-dimensional projective space and 
are i-spheres. Since they have the same universal cover x 3“^ 
and the same fundamental group Z 2 , they have the same homotopy 
groups but the Kiinneth formula implies that the Poincare polynomial 
of G = P^ X is Pg{x) = 1(1 + while the Poincare polynomial of 
H = P^ X is 1(1 -|- T^). Having different cohomology groups pre¬ 
vents the two graphs G, H to be homotopic. We see in this example, 
how useful it is to have a product in graph theory which shares the 
properties from the continuum. The point is that one does not have 
to reinvent the wheel in graph theory but that one can piggy-pack on 
known topology. 

From a practical point of view, the construction of the product only 
needs a few lines of code for a standard computer algebra system. The 
full computer code is given in detail later on. The product works for 
all hnite simple graphs: hrst construct the ring elements fniv) 

which are polynomial in xi,... ,Xn representing the vertices of G and 
1/1, • • •, I/m, the vertices of H, then take the product fafn and construct 
from this a polynomial in the variables xi,..., |/i,..., |/m, the new 

graph G X H hy connecting two monomial terms of a polynomial if 
one divides the other. For example, \f G = K 2 then fc = x + y + xy 
and if H = K 2 then fH = u + v + uv, we get /g/h = ux + vx + 
uvx + uxy + vxy + uvxy + uy + vy + uvy which encodes a wheel graph 
Ws with central vertex represented by uvxy. We see that the product 
naturally extends to chains, group elements in the ring and that it cor¬ 
responds to the product in the tensor product of the two rings. Every 
polynomial dehnes back a graph, but the later does in general not not 
agree with the original host graph. For example, take G = K 2 and the 
chain element ffx Pby + lf)xy, which dehnes in turn a graph which is 
the disjoint union of Ki (represented by 3a;) and K 2 represented by 
by -|- lf)xy. Chains are useful constructs because natural operations 
escape the class of graphs: examples are forming the boundary 5G or 
taking quotients G/A with a subgroup A of the automomorphism group 
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Figure 4. The product G of with (the two fac¬ 
tors are seen in the upper row) gives a 5-dimensional 
graph which has the same homotopy group than the 
product H of with (the two factors are seen in 
the lower row) which is also 5 dimensional. The graph G 
has 26 * 848 = 22048 vertices while H has 73 * 80 = 5840 

vertices. Both graphs have Z 2 as fundamental group and 
^2 

X as universal cover. But the cohomology of G 
is the one of 5^ and the cohomology of H is the one of 

by the Kiinneth formula. By the way, we constructed 
the projective spaces P^ by taking a rehned cross poly¬ 
top Ok which is a fc-dimensional sphere, then forming 
Ok X Ki which is large enough so that we can factor out 
the antipodal involution obtaining the projective space. 
This example illustrates the use of the product in con¬ 
crete graph theoretical constructions shadowing exactly 
the theory in the continuum. 


of G. Chains form a ring and have properties wanted from a geometric 
object: a dimension, a cohomology, a notion of homotopy, an Euler 
characteristic, a notion of curvature and these notions match the re¬ 
sults we expect from the familiar cases: Results like Gauss-Bonnet |16j . 
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McKean Singer and Hodge-De-Rham [19] , Poincare-Hopf UTl, Brower- 
Lefschetz m, Riemann-Roch m is a strong enongh theory so that 
it can be extended to higher dimensions), integrable geometric evo- 
Intions [231 E2], Riemann-Hurwitz m, Lusternik-Schnirelmann na, 
mirror the resnlts in the continnnm. See also [21] for an overview from 
the linear algebra point of view. These resnlts are more limited, if 
one considered graphs as one-dimensional simplicial complexes only, a 
common assnmption taken in the 20th centnry. As graph theory is a 
discrete theory, it can snrprise at first that the sitnation parallels the 
continnnm so well bnt there are conceptnal reasons for that, an ex¬ 
ample is non-standard analysis, an other is integral geometry. Some 
notions go over pretty smoothly, like cohomology or homotopy: others, 
like the notion of homeomorphisms for graphs needs more adaptations 

m- 


Abont the history: Kiinneth fonnd the formnla in 1921 HDl, where it 
was his dissertation nnder the gnidance of Heinrich Tietze. The paper 
pnblished in 1923 [32] is a linear algebra analysis which conld probably 
be simplihed considerably nsing Hodge theory. Kiinneth attribntes the 
Cartesian prodnct of manifolds to Steinitz (1908) (The definition is in¬ 
deed given in HU p.44 is probably the first appearance of the simplicial 
prodnct which onr dehnition is based on). The notion of manifolds was 
initiated by Poincare (1895), Weyl (1912), Veblen and Alexander [43] 
(1913) and Whitney [1^ (1936). It was Herbert Seifert who introdnced 
fibre bnndles in 1933 [22] • The Eilenberg-Zilber theorem of 1953 [40] . 
(Joseph Abraham Zilber was a Boston born Harvard gradnate (1943) 
and PhD (1963). Interestingly, the Eilenberg-Zilber theorem was an- 
thored 10 years before Zilber got his PhD degree nnder the gnidance 
of Andrew Gleason. More information abont the history of algebraic 
topology, see [2 EE] or the introdnction to [34] • The classical Carte¬ 
sian prodnct of graphs was introdnced by Whitehead and Rnssell in 
Principia Mathematica 1912 (it is [12] who spotted the constrnction 
on page 384 in Volnme 2 of that epic work). It was introdnced there in 
the context of logical relations, not so mnch graph theory. This histor¬ 
ical observation and many properties of the classical prodnct prodnct 
are discnssed in na. 

As far as we know, onr present paper is the hrst establishing a Kiinneth 
formnla for hnite simple graphs. There is a fnnctorial approach to 
Kiinneth for digraphs, where m use path cohomology to get a functor 
from digraphs to CW complexes, so that one can then use the contin¬ 
uum result for the CW complexes. Note however that in their case. 
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one gets Kiinneth only indirectly by constructing a CW-complex, take 
the product using the Cartesian embedding and then pulling the result 
again to graphs. In particular, the notion of dimension is also borrowed 
from the continuum as CW-complexes are topological spaces. A simi¬ 
lar thing could be done for geometric graphs, graphs for which the unit 
spheres S{x) are homotopy n-spheres. One can then build for each ball 
B (x) an open set in R"' and use this to build a manifold M. The open 
sets form then a nice good cover and the topology of the manifold is 
the same than the topology of the geometric graph. This construction 
is restricted however to geometric graphs. 

The missing Cartesian product bothered us for a while so that we 
decided to make a targeted search over several weeks (while procrasti¬ 
nating from an urgent programming job still in need to be hnished in 
geometric graph coloring), trying several possibilities and checking with 
the computer whether the construction works. The product described 
here was obtained by trying random things which look “beautiful”. 
The current product is attractive when seen algebraically because it 
becomes associative on the algebraic level. Until now, we would just 
take the usual product and then £11 out the chambers. When working 
with graphs, this is cumbersome, even when staying within the discrete, 
as it costs programming effort to build examples, like stellated higher 
dimensional cubes. The product which we propose here for graphs 
seems not to be known in the language of graphs. One could of course 
take the product |X| x |U| of the corresponding topological spaces as 
Kunneth did. After finding the product, we looked around whether it 
already exists: closest to what we do here appeared as an exercise in an 
algebraic topology course by Peter Tennant Johnstone at Cambridge 
[13]. More digging revealed that this simplicial product has been used 
by Eilenberg and Zilber in [iQ] (page 204), by De Rham [37| (page 191) 
and earlier by Steinitz in |TT] (page 44). But this product never made 
it to graph theory. One could also get a product by escaping to the 
Euclidean space. Our reluctance to use Euclidean stuff is not because 
we feel like Brouwer (who would even refuse to accept the infinity of 
natural numbers), but simply out of pragmatism: we want to build the 
structures fast on a computer without having to use Euclidean parts. 
Graphs are natural structures built into computer algebra languages 
and the Euclidean embeddings are not needed, when doing computa¬ 
tions; only if we want to see the graphs visualized or modeling tra¬ 
ditional geometric objects, the Euclidean embedding is helpful. Here 
is a self-contained full implementation of the new graph product in 
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“Mathematica”, a computer algebra system which uses graphs as fun¬ 
damental objects in its core language, structures which are dispatched 
from Euclidean embeddings, unless they are drawn. One procedures 
allow to translate a graph into a ring element and an other allows to get 
from a ring element back a graph. The graph prodnct is the product 
in the polynomial ring. 


Cliques [s_,k_] : = Module [{ n , t ,m, u , q , V=VertexList [ s ] ,W=EdgeList [ s ] , 1 } , 
n=Length[V]; ir^Length [W] ; u=Subsets [V, { k , k } ] ; q=Length[u]; 1={}; 
\^fcTable[{W[[j ,1]] ,W[[j ,2]]} ,{j ,m} ] ; If [ k==l, l=Table [ {V [ [ j ] ] } ,{j ,n}] , 

If [ k==2, l^V,Do[ t=Subgraph [s,u[[j]]]; If [Length [EdgeList [t]] = = 

Binomial [ k , 2 ] , l=Append [1 ,VertexList[t]]] , {j ,q}]]]; 1]; 

Ring [ s_ , a_] : = Module [{ v , n jUi, u ,X} , v=VertexList [ s ] ; n=Length [ v ] ; 
u=Table [Cliques [s,k],{k,n}] /. Table [ k—>a [[k]] ,{k,n}] ;n:^Length [ u ] ; 

X=Sum[Sum[ Product [ u [ [ k , 1 ,m] ] , 

{m, Length [u [[k,l]]]}] ,{1 , Length [ u [[ k ] ] ] }] , { k,m}] ] ; 

GR[ f _ ] : = Module [ { s = {} } ,Do[Do [ I f [Denominator [f [[k]] / f [[ !]]] = = ! SzSz k! = 1 , 
s=Append [ s , k—>1 ] ] , { k , Length [ f ] } ] , { 1 , Length [ f ] } ] ; 

UndirectedGraph [ Graph [ s ] ] ] ; 

GraphProduct [ si _ , s2 _ ] ; = Module [{ f,g,i ,fc ,tc}, 

fc=FromCharacterCode; tc=ToCharacterCode ; 

i [ 1 _ , n_ ]: = Table [fc[Join[tc[l], Int egerD igit s [k]+48]] ,{k,n}]; 
f=Ring [si , i [ ” a” , Length [VertexList[sl]]]]; 

g=Ring [ s2 , i [ ” b” , Length [VertexList [s2 ] ] ] ] ; GR[ Expand [ f *g ] ] ] ; 

NewGraph [ s _ ]: = GraphProduct [ s , CompleteGraph [ 1 ] ] ; 
example = GraphProduct [ CompleteGraph [ 3 ] , StarGraph [ 4 ] ] 

s._✓ 


The code of the above listing can be grabbed by looking at the source 
of the ArXiv submission of this text. 

Can we write the usual product in an algebraic way? Yes, if we have 
the ring element / and take as vertices the 0-dimensional simplices 
and as edges the 1-dimensional edges, we recover the old graph. If we 
take the product fg and take as vertices the pairs xy and as edges the 
triples xyz, we regain the standard Cartesian product. For example, 
iff = xy + X + y represents K 2 and g = uv + vw -|- wu + u + v + w 
represents K 3 , then the qnadratic and cnbic terms of fg are ux + vx + 
wx+uy + vy+wyuvx+uwx+vwx + uvy + uwy + vwy + uxy+vxy+wxy 
give the prodnct graph which has 6 vertices and 9 edges. 
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Figure 5. When invoking the product, like with 
GraphProduct[ConipleteGraph[3],StarGraph[8]] the pro¬ 
gram hrst computes the algebraic expressions in the ring, 
which are f = ai + a 2 + 0102 + 0,3 + OiOs + ( 12(^3 + 010203 
and g = bi b 2 & 1&2 “1“ ^3 + ^ 1^3 + ^4 + ^ 1^4 + ^5 + ^ 1^5 + 
bQ + bib^ + bj + bib-j + bs + bibg, containing 7 and 15 poly¬ 
nomial monoid entries, then multiplies them fg and de¬ 
rives from this again a graph Gfg = G x H. The product 
of the 2-dimensional triangle G and 1-dimensional star 
graph Ss is the graph G x H displayed in the hgure. It 
is 3-dimensional with vq = 105 vertices, vi = 203 edges, 
V 2 = 1182 triangles and ^3 = 1083 tetrahedra. Because 
both factors are homotopic to iFi, also the product is 
contractible and has the Poincare polynomial Pgxh{x) = 
1 . 
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2. Construction 

In this section, we describe that the vertices of a graph dehne a ring in 
which every element / can be seen as its own geometric object which 
carries cohomology, homotopy, Euler characteristic, curvature, dimen¬ 
sion and a Dirac operator. As we could run the wave or heat equation 
on such a ring element /, the chain / should be considered a geo¬ 
metric object with physical context. Unlike the category of graphs, 
the category of these chains has an algebraic ring structure, is closed 
under boundary formation as well as taking quotients which leads to 
orbifolds in the continuum. Other constructs in the continuum like dis¬ 
crete stratifolds are already implemented as suspensions of geometric 
graphs and are represented by classical graphs. We called them ex¬ 
otic in [28] , where we asked whether they might lead to exotic discrete 
spheres (a fc-dimensional graph with geometric unit spheres which is 
homeomorphic to a fc-dimensional homotopy sphere) and looked at dis¬ 
crete varieties in [2B]- Somehow, chains form are part of a list of class 
of structures which resemble structures in the continuum geometric 
graphs C stratifolds C varieties C graphs C orbifolds C chains. The 
product goes all the way to chains and can also be used to lift the 
classical Cartesian product to 1-dimensional chains: it is obtained by 
taking the product and disregarding higher dimensional parts of the 
chain: kind of projecting the result onto curves. 

We use multi-index notation ■ A hnite 

simple graph G with vertices Xi,..., dehnes a ring generated by all 
non-constant polynomial monoids x = Xk^^Xk^'' 'Xki and the chain ring. 
With a given orientation, graphs are always represented by functions 
of the form 'with Uj G {0,1} and 7^ 0- The star graph 

for example is x + y + z + w + xw + yw + zw. The ordering of the 
terms in the monoid allows in a convenient way dehne an orientation 
of the simplices in the graph G which means hxing a basis for discrete 
differential forms Q^{G). 

We could add constants to get a ring with 1 but we don’t yet see a use of 
the constants for now as we have no geometric interpretation for it. The 
Euler characteristic formula xif) = • • • 5 ~1) shows that the 

Euler characteristic of 1 is —1. The 1 element can’t be interpreted as 
the empty graph, because the empty graph is represented by 0, and is a 
(—l)-dimensional sphere with Euler characteristic 0. With a 1 element 
in the ring, one can produce terms like (l-|-a;)(l-|-j/) = 1 +x+y+xy, an 
object of Euler characteristic 0 or {2+xy){l+zw) = 2+xy+2zw+xyzw, 
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an object of Euler characteristic —6. 

To add an other footnote, we see that some chains represent already 
special differential forms, but these are forms taking values in the inte¬ 
gers. The triangle xyz + xy + yz + zx + x + y + z is a sum of a 2-form, 
a 1-form and 0-form which are all constant 1. In this discrete setup, 
geometric objects and differential forms are already very similar, as it is 
custom for quantum calculus setups, where Stokes is just the statement 
9 ) = (/) ^ 9 ) cind where the boundary operation is truly the adjoint 
of the exterior derivative and both geometric objects and forms are in 
the same function space. For classical differential forms and geometric 
objects, the geometric objects are distributions (as curves and surfaces 
for example are inhnitely thin) and differential forms are smooth or the 
dual setup used in geometric measure theory where differential forms 
are distributions and smooth functions are geometric objects. Its only 
on the level where both parts (geometric objects and differential forms) 
are represented by the same type of data that we have true symme¬ 
try, but then we are in a quantum setup. 

Every element / = dehnes a graph Gf: the vertices of Gf are 

the monomials of / and two monomials of / are connected by an edge 
if one is a factor of the other. The graph of / = Axy + 2x + y — 3x 
for example is the union of a line graph with three vertices where Axy 
represents the middle vertex, as well as a single Ki represented by 3a:. 
The graph G on the other hand dehnes the ring element fc = 
where x = Xj^... Xj^ is a simplex in G. The choice of the sign or 
permutation when writing down the polynomial monoid components 
corresponds to a choice of basis and is irrelevant for most considerations 
like for computing cohomology or when running discrete differential 
equations like [231122]. 

After hxing an orientation for each simplex, we get incidence matrices 
dfc. They implement the exterior derivative. They also determine the 
Dirac operator D = d + d*, and all such matrices are unitarily equiva¬ 
lent. Independent of the basis is the form Laplacian L = which has 
a block decomposition into Laplacians on /c-forms. The enhanced 
graph Gi = G/g is a rehnement of G and as we will see, if G is geo¬ 
metric, then the rehnement is geometric again of the same dimension. 
For general networks G we will show that the dimension of Gi can only 
increase, the reason being that higher dimensional parts will spawn off 
more “new vertices”. For a triangle G for example, the ring is Z[x, y, z] 
and fc = x + y + z + xy + yz + xz + xyz. This ring element dehnes the 
graph H = Gf^ for which the monomials x, y, z, xy, yz, xz, xyz are the 
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Figure 6 . The idea to associate a graph from a ring 
element can be done more generally. Let G be a graph. It 
dehnes a ring element / in the ring Z[xi,..., a;„], which is 
now abelian. We can form and look at the graph which 
this function defines. The hgures above show examples 
for the triangle for the circle G = G 15 and for the 
octahedron G = O. For G = in particular, we get 
the ring element f = a + b + c +ab + ac +be + abc and 
compute p = + 2a?b‘^c + a^b"^ + 2a%c^ + Aa^bc + 

2a?b + + 2a^c + P + 2a6^c^ + 4a6^c + 2ab‘^ + 4a6c^ + 

Qabc + 2ab + 2ac^ + 2ac + b‘^c^ + 2b‘^c + 6^ + 2bc^ + 2bc + 
which has has 23 monoidal summands. The incidence 
relations through factorization produces a graph with 23 
vertices. 

vertices. The divisor incidence condition leads to the wheel graph Wq 
which has the same topological features as G. Of course, the refine¬ 
ment process can be repeated, leading to larger and larger graphs with 
the same automorphism group. The sequence of graphs Gi, G 2 , G 3 ,... 
dehne so larger and larger rings generated by more and more variables. 
We will see that the dimension converges to an integer and that for 
large n, the graph is close to a geometric graph as the larger simplices 
will overtake all others. See Figure ( [IT| ). 

Given two graphs G, i/, dehne the rings Xc,Xh, take the product fg 
in the tensor product X oi rings and translate that product ring 
element back to a graph. The dimension of G is dehned inductively 
as 1 plus the average of the dimensions of the unit spheres and uses 
the foundation assumption that the empty graph has dimension — 1 . 
A graph is geometric, if every unit sphere is a homotopy sphere. A 
homotopy sphere is a geometric graph for which removing any vertex 
renders the graph contractible. Inductively, a graph G is called con¬ 
tractible, if there exists a vertex x such that S{x) and G\ {t} are both 
contractible, using the inductive assumption that the 1 -point graph is 
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contractible. Homotopy can be defined algebraically. Inductively, a 
ring element fa is contractible if there exists Xk such that both ring 
elements A{xi, ...Xk,.. .,Xn)) = f{x) - f{xi,... ,Xk-i,0,Xk+i, ...,Xn) 
and B{x) = f{xi ,..., Xk-i, 1, Xk+i, ■ ■ ■, Xn) are contractible. The Euler 
characteristic of G is x(G) = where Vk{G) is the num¬ 
ber of /c-dimensional simplices in G. More generally, the Euler 

characteristic of a ring element / is x{f) = • • • > ~1)- The 

Euler characteristic of the triangle G = x + y + z-\-xy + yz + zx + xyz 
for example is -(—3-1-3 — 1) = 1. The Euler characteristic of the chain 
?)X + by + 4xy is —3 — 5 -|- 4 = —4. The boundary of a ring element / is 
defined as Sf = ^xif, where dxi are the usual partial derivatives and 
the result is projected onto functions satisfying /(O) = 0. Due to the 
orientation assumption, this gives rise to sign changes. For example, 
5[xyz) = yz — xz + xy if the edges basis xy, yz, zx was chosen. We have 
66 f = 0. The exterior derivative on D is defined as dF{x) = F{ 6 x). 
As = 0, we have cohomology groups for any ring element. Let 
bk = dim(iL^(G)) be the Betti numbers. The Euler-Poincare formula 
x(G) = J2k=oi~^)^^k{G) follows from linear algebra. The boundary 
6 f{xiX 2 ...Xk) = ... ,xi,... ,Xk) is no more a graph in 

general. The star graph g = x + y + z + w + wx -|- wy -|- wz for ex¬ 
ample has the boundary 6 g = 3w — x — y — z which is a chain. The 
space of differential forms D is a direct sum D = ©D^, where is 
generated by functions supported on /c-simplices, polynomial monoid 
parts in the ring of degree k + 1. Since Q^{G) are finite dimensional, 
the maps dk are represented by finite matrices. They are called the 
incidence matrices and were considered by Poincare already for trian¬ 
gulations of manifolds. The cohomology groups H^{G) and more gen¬ 
erally H^{f) = ker((ifc)/im((ifc_i) are independent of the chosen signs, 
when defining the chain ring element fc. The matrix D = d + d* is 
called the Dirac matrix of /. Its square L = is the form-Laplacian. 
It decomposes into matrices Lk : D^(/) —t D^(/). By Hodge the¬ 
ory, the dimension of the kernel of Lk is the /c’th Betti number bk{G). 
Given a ring element / and a monoid part x in /, its unit sphere 
S{x) is the unit sphere of x in the graph defined by /. It consists 
of all monoids dividing x or which are multiples of x, without x. It 
is again a ring element. The dimension dimj(a;) is inductively de¬ 
fined as 1 plus the dimension of the unit sphere. The dimension of 
/ finally is the average of the dimensions of all monoids in /. The 
dimension of the triangle xyz + xy + yz + zx + x + y + z for exam¬ 
ple is 1 + [dim(S'(a;|/; 2 )) + dim(S'(a;|/)) + dim(S'(|/ 2 ;)) + dim(S'( 2 ;a;)) + 
dim(S'(a;))+dim(S'(|/))+ dim(S'(2:))]/7. We have for example S{xyz) = 
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xy + yz + zx + x + y + z which dehnes the graph Cq of dimension 1 and 
S{xy) = xy + x + y which is a line graph P 2 of dimension 1 etc. We see 
that the dimension of the triangle is 2. Of course, in the graph case, 
the dimension can be better computed on a graph level. The point is 
that the dimension extends to a nonnegative functional on the entire 
ring in such a way that the average of the dimensions of unit sphere 
S{x) of a monoids or generalized vertices of / is the dimension of / 
minus 1. 

We also need an inner product (/, g) on which is dehned as (a„x", 

= 'Yin ^nbn- It obviously satishes all properties of an inner product and 
especially dehnes a length |/| = \/(/, /). Of course, the incidence ma¬ 
trices d,d* are adjoint to each other with respect to this product. As 
graphs are special functions, we could use the inner product for exam¬ 
ple to dehne an angle between two graphs G,H on the same vertex set 
it is the arccos of the fraction Cov[G, H]/{a[G]a[H]), where Cov[X, Y] 
is the number of common simplices of G, H and cr[G] is the square root 
of the number of simplices in G. Lets summarize the main point; 

Proposition 1. a) Every graph G = {V,E) defines a ring element 
fc = Yx^y ® complete subgraphs of G. 

b) Every element f = Y^kX^ defines a graph Gj, where V are the 
monoid entries OkX^ in f, with \k\ >0 and where two entries are 
connected if one divides the other. 

c) If f comes from a graph G, then Gi = Gf is a graph for which 
the original simplices are the points and which has the same topological 
features than G and which additionally has a natural digraph structure. 

In other words, there is a functor f ^ Gf from the ring R to the cat¬ 
egory of directed graphs given by the division properties (even so we 
often forget about the directions) and that there is a second functor 
G —)■ /g from the category of undirected graphs to the ring. They 
are not inverses of each other but the cohomology agrees. For other 
functorial relations, see the recent paper [In- 

Remarks: 

1) Every ring element / also dehnes its own geometric object which 
has Euler characteristic, cohomology, dimension, homotopy as well as 
curvature. 

2) When forgetting about the anti-commutativity within the graph 
which is irrelevant for the graph product, the ring could be replaced 
with a more general integral domain. It would allow to see the graph 
product G X iL as an element in R[yi, ..., ym] where R = Z[xi,..., Xn\. 
This possibility can be useful when studying hbre bundles as one can 
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work in a ring of the fibre graph. 

3) The fact that the product of two graphs is obtained by writing the 
graphs algebraically using different generators and producing from it 
again a graph: 

G X H = Gf^xfH 

is not unfamiliar to us. If we take the Cartesian product of two spaces, 
we use different variables for the different directions. If we don’t take 
new variables and take / ■ / in the ring, then this is in general a chain. 
For a triangle G = xy + yz + zx + x + y + z + xyz for example, we get 
(using x^ = y"^ = z'^ = 0) the chain G ■ G = 2xy + 2yz + 2zx + 6 xyz. 
The graph which belongs to this chain is the star graph S 3 as 6 xyz is 
the central vertex and the others the outer points (as they divide the 
central point). 

The imposed Pauli principle imposed by anticommutativity xy = —yx 
is irrelevant for the Euler characteristic, both for chains as well as for 
the graphs Gf derived from the ring elements /. 

4) The graph G x K 2 can be seen as a self-cobordism of Gi with itself 
as it is a graph of one dimension more which has two copies of Gi 
as boundary. It is in general true that any geometric graph is self 
cobordant to itself as we can sandwich two copies G with a completed 
dual graph G. but for Gi we don’t have to work on a construction. It 
is given. 

5) Denote by 11^^ the projection onto the linear subspace generated by k 
until /-dimensional simplices. As pointed out before, we can recover G 
from fc by building h = Hoifc and then building Gh which is G. With 
h = 1123(/g/h), then Gh is the classical standard Cartesian product of 
G and H. 


3. Some geometry 

The results for the graph product mirror results in the continuum. 
First we look at some basic constructions which deal with the notion 
of homotopy sphere or simply fc-sphere in graph theory. The dehnition 
of a /c-sphere in graph theory is recursive: a /c-sphere is a /c-dimensional 
geometric graph for which every unit sphere S{x) is a (/c — l)-sphere 
and such that after removing any of its vertices, we get a graph which is 
contractible. A k-haAl is a /c-dimensional geometric contractible graph 
with boundary which has a (/c — l)-sphere as its boundary. We call 
the interior of a ball the part of B which is not in the boundary. A 
suspension SG of a graph G is the join G x Sq, a double pyramid 
construction: add two new points x, y and connect the points to all 
the vertices of G. A pyramid construction itself is the join G -k Ki. 
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From the Cartesian product, we have construct joins by just building a 
product and identifying some variables in the algebraic representation 
of the graph. Examples are given at the end. 

Lemma 2 (Suspension). The join G'kKi of a k-sphere G is a {k + 1)- 
ball. The suspension G -k K 2 of a k-sphere G with a 0-sphere K 2 is a 
{k + l)-sphere. 

Proof, a) By dehnition, the boundary of Sk Ki is S which is a sphere. 
Also, the graph S k Ki is contractible. 

b) Removing the second point y produces the ball S k Ki hj a). □ 

Remark: 

1) More generally, as in the continuum, and shown below, the join 
k of two spheres is a sphere 

2) Also as in the continuum, the dehnition of the join needs the prod¬ 
uct as the join Ak B is a quotient of A x B x I. It generalizes that 
the 3-sphere can be written as k which has an interpretation of 
gluing two solid tori along a torus. 

Examples: 

1) The wheel graph Wn is a 2-ball. It is the join Gnk Ki, where Gn is 
the cyclic graph with n vertices. 

2) The /c-dimensional cross polytope is S 2 kS 2 k ■ ■ -kS 2 , where we have 
{k -|- 1) factors. The square C 4 is equal to S 2 k S 2 , the octahedron is 
S 2 k S 2 k S 2 etc. 

Lemma 3 (Glueing ball). Assume Bi are k-balls with boundaries Si 
and assume that R 3 = Ri fl B 2 is a {k — 1 ) ball with boundary S 3 = 
SiA S 2 . Then B = BiU B 2 is a k-ball with boundary S (Z Si U S 2 . 

Proof. This is proven by induction with respect to k. There are four 
things to show: 

a) B is contractible. 

b) every unit sphere in the interior of R is a (/c — l)-sphere. 

c) every unit sphere in S' is a (/c — 2)-sphere. 

d) when removing a vertex from S, we get a contractible graph. 

For a) take a point y in B 3 . We can retract everything in Bi,B 2 , B 3 to 

y- 

For b) we only have to look at a vertex x in B 3 . The unit ball B{x) 
decomposes For c), we only have to look at a vertex 2 ; in S 3 and see 
whether its unit sphere in S' is a homotopy sphere. For d), we can 
retract a pointed part to S'3. □ 
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The next statement is the discrete analogue of the classical statement 
that the boundary of the product of two balls d{Bi x B 2 ) is a sphere 
and that it can be written as as dBi x B 2 + B 1 X dB 2 which is the union 
of two solid tori glued at a torus. Also in the discrete, we can use the 
intuition from the continuum: 

Lemma 4 (Cylinder lemma). If Bi are ki-balls with {ki — l)-spheres 
Si = 6 Bi as boundary, then (i?i x S' 2 )U(i ?2 x S'!) is a {ki + k 2 — l)-sphere 
provided Bi x 82 ^^ B 2 x Si d Si x S 2 - 

Proof. Use induction with respect to dimension. The union is a graph 
of dimension ki + k 2 — 1. A unit sphere S{{x,y)) is in Bi{x) x S 2 {y) 
or Si{x) X B 2 {y) whose intersection is Si{x) x 5*2(i/). □ 



Figure 7. The Cylinder Lemma tells that “the union 
of complementary solid tori glued together at a torus is 
a sphere.” To the left, we see Bi x Sq U Sq x Bi which 
consist of the union of a bottom-top and left-right pair. 

The intersection Bi x SqDBi x Sq is equal to Sq x Sq. The 
second example seen to the right shows B 2 X SqUSiX Bi 
which is the union of top-bottom discs and a mantle. The 
intersection B2 x SqP Si x Bi agrees with Sq x Si which 
is the union of two circles. We are interested in these 
type of spheres because unit spheres of product graphs 
are of this “Hopf hbration” type. 

Examples: 

1) If Bi are two line graphs with boundary Si, then Bi x 5*2 is a union of 
two line graphs. Similarly B2 x Si is the union of two line graphs. The 
union is a square, the intersection consists of the four points Si x Si. 

2) If Bi is a ball with 2-dimensional sphere Si and B 2 is a 1-dimensional 
ball with 1-dimensional sphere S2, then Bi x S2 is the mantle of the 
cylinder and B2 x Si are the top and bottom cover. 

3) If Bi is 1-dimensional and Bq is 0-dimensional, then dBi x Bq is a. 
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0-sphere. 

4) If Bi, B 2 are two-dimensional balls, we can see dBi x B 2 + B 1 X dB 2 + 
dBi X dB 2 is the union of two solid tori glued along a 2-torus. This Hopf 
fibration is classically given as the split of = {{z,w) G \ \z\‘^ -|- 
[tcp = 1} into two solid tori Bi = {|; 2 | > |tc|} and B 2 = {|;2| < |tc|} 
intersecting in the 2-torus {| 2 :| = |tc| = 1}. 

4. De Rham cohomology for graphs 

Among various other flavors of cohomologies, there are three equiva¬ 
lent cohomologies for compact n-manifolds: simplicial cohomology, de 
Rham cohomology and Cech cohomology. For simplicial cohomology, 
the manifold is triangulated into hnitely many n-simplices leading to a 
differential complex. In de Rham cohomology, one works with the com¬ 
plex of differential forms, partial derivatives tap into the local product 
structure of the manifold, for Cech cohomology, the manifold is cov¬ 
ered with a hnite cover of open sets so that the nerve graph determines 
the cohomology. Each of the these cohomologies have advantages over 
the others: simplicial cohomology is the computer science or combi¬ 
natorial point of view which sees space as a mesh of small simplicial 
building blocks, the de Rham cohomology is the analysis or calculus 
approach, which taps into the bag of techniques used in calculus. In 
this flavour, the basic building blocks are cubes obtained from a local 
product structure and the different directions are accessed with partial 
derivatives which we think of the exterior derivatives in each factor. 
The Cech cohomology hnally is the homotopy or topologal point of 
view which relies on the fact that cohomology is more robust and tran¬ 
scends dimension. Simplicial cohomology is straightforward and the 
simplest. De Rham cohomology requires local charts which are prod¬ 
ucts and taps into the differential structure of the manifold, allowing 
for an efficient computation of cohomology. Cech cohomology hnally 
is a hexible variant which illustrates best the homotopy invariance of 
cohomology. It also can lead to considerable complexity reduction as 
on can for example retract a space to a much smaller dimensional set. 
A solid torus for example can be retracted to a circle. The equivalence 
of simplicial cohomology with de Rham cohomology is due to de Rham. 
Later proofs of this theorem use the equivalence of Cech cohomology 
with simplicial cohomology. 

All three cohomologies have analogue constructions in graph theory. 
It is our goal to introduce the analogue of de Rham cohomology and 
give a discrete analogue of the de Rham theorem. Of course de Rham 
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cohomology emerges only in the discrete if one has a product which is 
compatible. 

The simplicial cohomology of a graph G is the clique cohomology of G 
using the Whitney complex of all complete subgraphs. It is the oldest 
and has already been considered by Poincare, even so not in the lan¬ 
guage of graphs. 

The Cech cohology is the cohomology of the nerve graph of a good 
cover, where “good” in the discrete means that the nerve is homotopic 
to G. Cech cohomology has hrst been considered for graphs in izg, 
a paper which proposes a notion of what a continuous map between 
graphs is, which is more tricky than one might think at hrst, as clas¬ 
sical topology badly fails as the topology generated by the distance is 
discrete making it unsuitable. What is important for a good notion of 
continuity is to merge homotopy with dimension. The equivalence of 
discrete Cech cohomology with discrete simplicial cohomology is there 
by dehnition just relying on the fact that homology is a homotopy in¬ 
variant. 

What was missing so far in graph theory is an analogue of de Rham co¬ 
homology, where “cubes” rather than “simplices” play the fundamental 
role. But one can not really look at a de Rham cohomology for gen¬ 
eral graphs, if one does not have a Cartesian product for which there 
is compatibility. A de Rham type cohomology has been mentioned in 
m, where a generalized path cohomology introduced by [S] is consid¬ 
ered for digraphs. As their approach uses a functor from graphs to CW 
complexes and pulls back results from the product of CW complexes 
to digraphs, there is no relation with what we do here. Other takes 
on discrete de Rham cohomology study discrete notions for numerical 
purposes |12] or |2]. Our approach to de Rham cohomology is purely 
combinatorial and restricted to hnite constructions. 

As in the continuum, the de Rham complex for a product graph G x 
H can use some derivatives also in the discrete, but this is merely 
language: while in simplicial cohomology, we write 5 {xiX2 ■ ■ ■ Xn) = 
^^(— 1 )^X 1 ... hj... Xn for the boundary of a simplex X 1 X 2 ... x„ in the 
de Rham approach, we can write 5x = Yli 9xiX in the algebra which 
is the same thing. It just uses the derivative notion in a formal way. 
The de Rham connection will be needed in the Kiinneth connection, 
where we look at the kernel of form Laplacians. Kiinneth will then be 
quite obvious. Without linking de Rham with simplicial cohomology 
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the relation is nontrivial, as the dimension of the space of differential 
forms on the product graph is much larger than the product of the 
dimensions of the space of differential forms on the factors. Already 
Kiinneth had to work though such difficulties and needed dozens of 
pages of linear algebra reductions to tackle the issue. Our situation is 
also different in that we look at the product of two arbitrary networks 
G, H which by no means have to be geometric. 

The analytic de Rham approach allows the derivative d{fg) to be writ¬ 
ten as a sum of products {df )g + {—iy^^f{dg) which is Leibniz formula 
and which reduces the exterior derivative of the product to the exte¬ 
rior derivative of the factors in the same way than the gradient, curl or 
divergence reduces the exterior derivative to partial derivatives, which 
are the exterior derivatives in the 1-dimensional factors. To illustrate 
this with school calculus: inhnitesimally, the curl of F = [P, Q) is 
Qx ~ Py As a graph theorist we look at P (for hxed y) as the 1-form 
restricted to the first coordinate x which means a function on edges 
of the first graph and at Q (for hxed x) as the 1-form restricted to 
the second coordinate which corresponds to edges in the second graph. 
The curl is so a line integral along a square. When reducing this to 
simplicial cohomology, the square needs to be broken up into triangles. 
Our product does that very explicitly even so it needs some care as the 
tensor product of hnite dimensional algebras has a completely different 
dimension in general than the Cartesian product. Kiinnneth needed 
dozen of pages of rather messy linear algebra reductions to achieve 
this. The language of polynomials and the de Rham connection will 
allow us to make this more clear. 

When we start with a triangulated picture, there are no more two dis¬ 
tinguished directions present. In graph theory, we only can distinguish 
two directions if we look at the product G x H of two graphs. As for 
manifolds, this structure could be allowed to be present locally only; 
its important however that a product structure must be present before 
we can even talk about discrete de Rham cohomology. As mentioned 
before, there is the possibility to see a graph as a triangularization of a 
manifold and use the Euclidean product structure to emulate a discrete 
de Rham cohomology. Notice however that this does not tickle down 
to the discretization. The relation would only exist functorially and is 
pretty useless when working with concrete networks. We will not leave 
the discrete realm and show that /c-forms on the product space can be 
related to products of forms in the two factors. And also, we do not 
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only work with geometric graphs, which can be seen as discretizations 
of manifolds; we work with general hnite simple graphs. 

A hnite simple graph G naturally comes with a simplicial complex, 
given by the set of all the complete subgraphs of G. This so called 
Whitney complex can be encoded algebraically in the ring of polyno¬ 
mials. For a triangular graph for example, we have the ring element 
g = xyz+xy + yz + zx + x + y + z, where the choice of the orientations of 
all the simplices is done arbitrarily. The ring element g in turn dehnes 
a new graph Gi in which the polynomial monoids form the vertices and 
two vertices are connected if one divides the other. It is important that 
it actually can be seen as a digraph, the direction is given which part is 
a factor of the other. In some sense, going from G to Gi frees us from 
having to chose an arbitrary orientation as the structure is now built 
in. By the way, Gi appears to have other nice features like having the 
Eulerian property allowing therefore a geodesic dynamical system |28j . 

In our triangular graph example G, we get a graph Gi with 7 vertices 
because there are 7 complete subgraphs of the triangle. The graph Gi 
is a wheel graph which shares the topological and cohomological prop¬ 
erties with the triangle. It is even homeomorphic to the triangle in the 
sense of [2^ as we can End a 2-dimensional open cover whose nerve is 
the triangle. Indeed: 

the Cecil cohomology of Gi is equivalent to the graph cohomology of 

G. 

The observation that G and Gi are homotopic proves that the Cecil 
cohomology of a graph with respect to a good cover in the sense of 
IZH is the same than the graph cohomology. As in the continuum, the 
relation between simplicial and de Rham cohomology is not completely 
obvious as we will just see. 

Now, if we take the product G x if of two graphs, like for example 
two complete graphs G = K 2 ,H = K 2 (illustrated in Figure ([^), 
where the ring elements g = {x + y + xy) and h = u + v + uv en¬ 
code the graph, then the product G x H is encoded by a ring element 
gh = XU + XV + xuv + yu + yv + yuv + xyu + xyv + xyuv which 
by looking at division properties of the polynomial monoids produces 
the graph G x H with 9 elements. It is the wheel graph J = Wio, 
which is a discrete square with 16 edges and 8 triangles. The Lapla- 
cian L = {d + d*y with the usual exterior derivative is a block matrix 
decomposing into a 9 x 9 block Lq, a 16 x 16 block Li and a 8 x 8 
block L 2 . The algebraic representation of f{J) is a polynomial element 
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with 9 + 16 + 8 = 33 monoid terms. The graph Ji associated with 
/(J) would have 33 vertices already. From topological, algebraic or 
homotopical considerations, the graphs J and Ji are equivalent: they 
are homeomorphic, they are homotopic and have the same cohomology 
and dimension. 

A boundary on the product is dehned by the Leibniz formula: 
SUg)^(6f)g + {-lf'f(6g) 

which then defines an exterior derivative d. As for Cech cohomology, 
we will see that there is an advantage in that we can work with a 
smaller dimensional vector spaces: the spaces Q^{G) and Q^{H) and 
even the tensor product Q’^{G) ® Q^{H) have in general smaller di¬ 
mension than the vector spaces fl^{G x H). See Figure Q. Having 
equivalence of cohomology can be a blessing when doing computations. 

This will be useful when looking at discrete manifolds: De Rham coho¬ 
mology can be used more generally also when gluing product graphs to 
build hbre bundle, this product structure does not have to be global. 
It is the same situation as for manifolds have in general only locally 
neighborhoods which can be written as products. 


Theorem 5 (Discrete De-Rham Theorem). The de Rham cohomology 
for a product graph G x H with boundary operation defined by the Leib¬ 
niz rule d{fg) = {df)g + {—iy-^^f{dg) using the exterior derivatives d 
on each factor G and H, is isomorphic to the cohomology on G x H 
defined by the exterior derivative d on the graph G x H. In short, 

H*,j,{G xH)^H*{GxH). 

After de Rham [3^, new proofs of the de Rham theorem in the con¬ 
tinuum were given by A. Weil (1952) [H] and H. Whitney |16]. Weil 
had outlined a proof already in 1947 in a letter to H. Cartan, a letter 
which initiated Cartan’s theory of sheaves. Weil’s argument used a 
staircase argument in a double complex. The textbook proof of Bre- 
don [35] is worked out in more detail in [35]. See also [36113119]. In 
the discrete, we can directly construct the chain maps 0 and 0. As 
we will see below, our goal will be to find a one-to-one correspondence 
between harmonic forms on the large simplex Laplacian L = (d + d*)^, 
which is a v{G x H) x v{G x H) matrix and the de Rham Laplacian 
L = {ddR + d'^j^Y, which is a (u(G) •u(id)) x (u(G) •u(id)) matrix. This 
is done by a rather explicit chain homotopy. 
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Figure 8. The product oi G = K 2 = x + y + xy with 
H = K 2 = u + V + uv is the wheel graph G x H = 

XU + XV + xuv + yu + yv + yuv + xyu + xyv + xyuv. It 
has a natural digraph structure in which the direction is 
given by which part divides which. A form / in hi* [G) is 
a function f = ax + by + kxy, where (a, b) represents the 
scalar function, the 0-form and k is the value of the 1- 
form. Similarly G Q*{H) is a function = du+ev+luv. 

The product fg G hi* (G) ® fl* (H) gives a scalar function 
on the 9 vertices of G x if. It dehnes an element F = 

^J^fg) in r2*(G X H) as the digraph structure dehnes an 
orientation on the simplices a and F{a) is the sum of 
the values of fg on the vertices. The reverse chain map 
0 is obtained by hrst dehning from F a function h on 
the vertices of G x if, then hnding the functions / and g 
by solving a system of equations. Now 4’i^) = fd is an 
element in Q* (G) {H). We see that differential forms 

and geometric objects are treated similarly, as is typical 
in quantum calculus setups and that the Eilenberg-Zilber 
becomes combinatorial and explicit. 

Proof. Assume G, H have vertices Xi,..., and i/i,..., y^. The de 

Rham cohomology on the ring generated by functions fg is dehned 
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by using the exterior derivatives on each product using the Leibniz 
rule. We have to relate it with the Whitney chain complex x H) 
which features functions on the simplices oi G x H. The linear spaces 
X H) and Vt^{G x H) have different dimensions VkiG x H) and 
Ylim+i=k '*^m{G)vi{H) so that (j) is not the inverse of -ip (they are chain ho¬ 
motopic only, which essentially means that they are equivalent modulo 
coboundaries on each side). In both cases, denote by d : —)■ 

and by d : hi"' —)■ the exterior derivatives. We will construct two 

linear maps 0, ip and check that they are chain maps: dip = ipd and 
d(p = (pd. This will establish the isomorphism of cohomology H*{Gx H) 
to X H). 

Construction of ip : —)■ 

Start with an /c-form F G be written as f{x)g{y) in the 

variables Xi,..., Xn, yi,..., ym, h dehnes a function on the vertex set 
of G X TT. Given a /c-simplex a in the graph G x H. Assigned to it the 
value of the vertex which belongs to a monomial of degree k. If there 
is one, it is unique. If there is none, assign the value 0. 

Construction of 0 : —)■ 

For every /c-form F in Q^{G x H) we build a polynomial (p{F) = 
f{x)g{y) as follows: the function F dehnes a value F to the vertices of 
G X H hj averaging the F values of the /c-simplices hitting x. We have 
now a function F on the vertices of G x H, written as a polynomial 
G{x, y) in the variables Xi,... ^x^yi, ■ ■ ■ ,ym- One can now add a poly¬ 
nomial P of smaller degree so that G{x, y) + P{x, y) can be factored as 
f{.x)g{y). 

0 is a chain map: ipd = dip: 

Proof: Given a. {k — l)-de Rham form F = f{x)g{y), build the k- 
form dF using the Leibniz rule, then use ip to get from this a /c-form 
ipdF on G X H. By linearity, we can assume that F = Ax^y^ with 
|i| -|- |j| = k — 1. The function dF{x,y) takes the value ±A on each 
polynomial y^ with \i'\ + \j'\ = k. Now, ipdF is a /c-form on G x H 
which takes the valne FA on those vertices y^ . Bnt ipF is a /c — 1 
form which assigs the value A to every simplex having x^y^ as a vertex. 
Now dip{F){x, y) is exactly A if xy is of the form y^ . 

0 is a chain map: (pd = d(p: 

Proof: given a k — 1-form F on the product graph G x K, we build 
the {k + l)-form dF and produce from this a polynomial ip{dF) in the 
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de Rham complex. We check that this is the same than dip{F). by 
linearity we can look at a /c — 1-form F which assigns the value A to 
a (/c — 1) simplex a of G x H. Now take a /c-simplex which has a in 
its boundary, then dF{a) = A. But now, (^{dF) assigns the value A 
to the vertex 2 ; which is in a and has degree k. Now, 4>{F) assigns the 
value A to the vertex z' which is in a and has degree k — 1. This is 
connected to so that d(l){F) is the same value. □ 

It would be nice to have a more intuitive understanding of the function 
6 : Q^j^{GxH) —)■ Q’^~^{GxH) satisfying the chain homotopy condition 

(ptl) — ip(f) = dO + 9d . 

The map 9 is related to the lower degree polynomial P added to F 
so that F{x,y) + P{x,y) = f{x)g{y). This involves solving a system 
of linear equations. The dimension of the space of scalar functions on 
Gx H is v{G) X v{H), where v{G) is the number of simplices in G. The 
space of functions f{x)g{y) has dimension v{G)+v{H). The image of 9 
consists of all lower degree polynomials (as usual satisfying /(0,0) = 0) 
which is a space of dimension n = (n(G) — l)v{H) + v{G){v{H) — 1). 
To solve for F{x, y) = f{x)g{y) we have to solve a system of equations 
for n variables and this solution gives us P. 


n 


k-l 

dR 



n 


I 

fc-i 



Ir 


Examples: 

1) If G is a forest with n trees and H is a. forest with m trees, then 
the product graph G x H has n*m components. The Betti vectors are 
b{G) = (n, 0,...), b{H) = (m, 0,...) and b{G x H) = {nm, 0,0,...). 
Every harmonic form on the product space G x H is of the form 
f{x)g{y), where f,g are both locally constant. Because the number 
of vertices in G x if and 0-simplices in G x if are the same, we don’t 
have to translate between simplicial and de Rham 0-forms. 

2 ) Let G = H = K 2 , then G x if is a 2-dimensional ball with 9 ver¬ 
tices {xu,xv,yu,yv,xyu,xyv,xyuv,xuv,yuv}. Given a 1-form F in 
the product graph G x H. It is a function on the edges of G x if. 
Given a vertex which is the product of a vertex and edge like xuv, we 
assign the sum over all edge values hitting the vertex. 


Remarks: 

1) In the complement of the set R of chains which are a product of 
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Figure 9. An example illustrating the chain homotopy 
in the case G = K 2 , H = K 2 '. a 2-form in x H) is 

given as a function f{xy)g{uv) as there are no 2-forms in 
G, H. The chain map "0 constructs from this a 2-form in 
X H), that is a function on the triangles of G x H, 
as follows: the function F{x,y,u,v) = f{x,y)g{u,v) = 
I2xyuv is a function on the vertices of G x if which is 
12 on the center point and zero else. The 2-form ‘ip{F) 
assigns now the value 12 to all triangles having xyuv as a 
vertex. The backwards chain map 0 assigns to the node 
xyuv the average of the values of the 2-form on adjacent 
triangles. There is no need here to add lower degree 
polynomial dO as this already factors f{x)g{y). 


graphs, the equivalence does hold any more. Take the chain / = 7xy + 
3x -|- 5y. Its cohomology is different from the corresponding graph Gf 
which is the 3 vertex graph without edges. 

2 ) The graph G/ of an abstract chain / has as many vertices as there 
are simplices in /. The simplicial complex of Gj therefore is already 
large. But the linear space of all fc-forms on G x ii is huge: if v{G) is 
the number of simplices of G and v{H) the number of simplices in H, 
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G 

Figure 10. An example illustrating the chain maps in 
the case G = K 2 , H = K 2 : a 1-form F in x H) = 

(g) © r2^(G') © is represented alge¬ 

braically as F = cxy{au + bv) + Guv{Ax + By) = 
{5uv){3x) = 15xuv. By looking at the coefficients, we 
get a function on the vertices of G x if. We have now 
assigned values to the four “1-dimensional vertices” rep¬ 
resented by monomials of degree 3, in this example to 
xuv only. The 1-form ^|JF on G x if assigns values to all 
edges emanating from xuv and assigns there the value 15 
attached to the node yuv. In this case we already have a 
factorization / = 3a:, g = 5uv. In general, we have to add 
a gradient to d6 in such a way that the chain homotopy 
— %l)(f) = dO + 9d holds. 

then Binomial(n(G) • v{H), k) is the dimension of the set of fc-forms on 
G X H. The de Rham complex is more manageable on a computer. 

















































THE KUNNETH FORMULA FOR GRAPHS 


33 


5. Results 

Theorem 6 (Geometry of product). If G,H are geometric graphs of 
dimension k, I then G x H is a geometric graph of dimension k + 1. 

Proof. We have only to show that S{{x,y)) is a homotopy sphere. De¬ 
pending on {x,y) the graph is the union of two cylinders intersecting 
in a lower dimensional graph which is part of the sphere S{x) x S{y). 
Use the cylinder lemma. □ 

Corollary 7. If G, H are geometric spheres of dimension k, m then 
the join G x H is a geometric sphere of dimension n = k + m + 1. 

Proof. It is a graph of dimension n -|- A; -|-1 as it is a quotient of a graph 
G X H X Ln, where Li is a line graph with I > 2. This assures that 
for all points which are not affected by the identification, we have a 
geometric unit sphere of dimension n — 1. For the vertices y which 
were subject of identihcation, the unit sphere Siy) is a union of two 
(n — l)-dimensional cylinders which by the cylinder lemma is a sphere. 
Besides verifying that every vertex y has a unit sphere S{y) which is 
a (n — l)-sphere, we have also to check that taking away a vertex x 
from Gx H, leads to a remaining graph which is a contractible ball of 
dimension n which has a (n — l)-dimensional boundary. □ 

Take Gi = G x Ki. Given a vertex t G Gi. It corresponds to a simplex 
in G. The set S{x) consists of all simplices in G containing x united 
with the set of all simplices inside x. 

Examples: 

1) Let G be a cyclic graph Gn- Lets look at the sphere of the vertex 
X which is now also a vertex in Gi. It becomes now the set of two 
vertices {yx, xz} in Gi corresponding to two edges in G. This is a 0 
dimensional sphere. 

2) Let G = Gn again and let e = xy he an edge in G which becomes 
now a simplex in Gi. Its sphere consists of the two simplices {x,y} 
which correspond to two vertices in G. The sphere is S'o, again a geo¬ 
metric graph. 

3) In the case d = 2 like if G is an ecosahedron, assume that S{x) = 
{yi,y 2 , ■ ■ ■ ,yk) is the old unit sphere of a vertex x. The new sphere 
of X in Gi contains the points xyk which correspond to edges in the 
old graph G. The triangles xyuyk+i are the new connections in that 
sphere. We have again a geometric graph. 4) Now assume that e = xy 
is an edge in the old two dimensional graph G which becomes a vertex 
in Gi. The sphere S{e) consists of the vertices t, y which are the third 
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points in the triangles xyzk as well as the two triangles xyzk- These 
two original vertices and triangles form 4 points in the new graph Gi 
which form a circular graph This is our unit sphere. 

5) Now assume that t = xyz is a triangle in the two dimensional graph 
G. It becomes a vertex in the new graph G\. The unit sphere S{t) 
consists of all simplices x, y, z, xy, yz, xz which is a 6-gon. Examples 
3)-5) together show that that Gi is always geometric if G is two di¬ 
mensional. The unit spheres always are G4 and Gq. In particular, Gi 
is Eulerian. 

6 ) In the case d = 3, if a; is an original vertex, then the sphere consists 
of all edges and triangles in the original sphere by induction this is 
S{x) X Ki. If e = {x,y) is an edge. Then the sphere consist of the 
points X, y as well as all triangles and tetrahedra containing e. This is 
a double suspension of a 1-dimensional sphere and so 2-dimensional. 
If f = (a, b, c) is a triangle, then the unit sphere is a suspension of the 
6-gon a, b, c, ab, ac, be. 


A theorem from 1964 assures that only spheres are manifolds which are 
joins [33]. Since geometric graphs naturally dehne compact manifolds, 
where unit balls are hlled up to become the charts, this holds also for 
geometric graphs. Intuitively, the reason is clear. In order that the 
unit spheres of the identihed end points are spheres, we better need 
the factors A, B of the join A-k B to be spheres, so that G is of the 
form A-k B where A, B are spheres. 


Corollary 8 (Kwun-Raymond). If G is a n-dimensional geometric 
graph which is the join A-kB for two other graphs, then G is a homotopy 
sphere. 

Now, we come to the main result as announced in the title: 

Theorem 9 (Kiinneth). Given two finite simple graphs G,H. The 
cohomology groups of G x H are related to the cohomology groups of G 
andH by 

H\G xH)= H\G) ® W{H) . 

i-\-j=k 

Proof. We know from Hodge theory that the de Rham Laplacian L = 
{d -|- d*Y = dd* -|- d*d of G restricted to fc-forms has a kernel ker(Lfc) 
of dimension is H^{G). 

Weha.ved{fg) = {df)g+{-iyi^ f{dg) and d*d{fg) = {d*df)g+f{d*dg) + 
{df){d*g) {d*f){dg), so that 

LUg) = {Lf)g + f{Lg) + 2{df){d*g) + 2{d*f){dg) . 
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Hodge theory also assures that that Lf = 0 is equivalent to df = d* f = 
0, so that if /, g are harmonic forms in G or hf respectively, then fg is a 
harmonic form in Gx H. We see that we can construct the cohomology 
classes of the product W{G) ® H^{H) from the cohomology classes of 
the factors G, H. It is obvious therefore that H^j^{G x H) contains the 
vector space ® 

Now assume that fg is harmonic in G x H. We want to show that 
Lf = 0 OT Lg = 0. Again we use that the kernel of L is the intersection 
of the kernel of d and the kernel of d*. Now look at the equation 

L{fg) = {Lf)g + f{Lg) + 2{df){d*g) + 2{d*f){dg) = 0 , 

where / is a i-form and g is a j-form. When we look at the terms, then 
/, {Lf) are polynomials in Xk of degree i + 1 and Lg, g are polynomials 
in Hk of degree j + 1. But {df){d*g) consists of polynomials of degree 
i + 2 in Xk and of degree j in yk and {d* fg) consists of polynomials of 
degree i + 2 in the Xk and of degree j + I in yk- Because the degrees 
are different, we see that 2{df){d*g) = 2{d*f){dg) = 0. Now, either 
df = d*f = 0 which implies Lf = 0. But then, f{Lg) = 0 so that also 
{Lg) = 0. An other possibility is dg = d*g = 0 which implies Lg = t] 
and in turn implies Lf = 0. A third possibility is df = t) = dg = t]. 
But from Lfg + fLg = 0 follows then {dd*f)g + fdd*g = 0. Now take 
the inner product with fg to get {f, dd*f) {g, g) + {f, f) {g, dd*g) = 0 
which implies that f = g = 0, demonstrating that the third case is not 
possible. We have seen that Lf = 0 and Lg = 0 and verihed that every 
kernel element in L{fg) necessarily has to have the property that both 
Lf = 0 and Lg = 0. 

Having established that the de Rham cohomology oi G x H satishes 
the Kiinneth formula and that the de Rham theorem assures that it is 
equivalent to graph cohomology, we are done. □ 

Remark: In the continuum, the equivalence of the chain complex 
xR) with the tensor product on which the exterior 

derivative is given by the Leibniz formula is subject of the Eilenberg- 
Zilber theorem from 1953 [10]. They use the nerve of product coverings 
and so Cech type ideas to verify the equivalence of the cohomologies 
and establish so the Kiinneth formula more elegantly. The above com¬ 
putation just reflects the fact that the tensor product of two chain com¬ 
plexes with exterior derivative given by the Leibniz formula is again a 
chain complex and that we have an Eilenberg-Zilber formulation as in 
the continuum, where two chain complexes A, B are chain homotopic, 
if there are chain maps 0 : A —>■ B,'tpB —)■ A such that the composi¬ 
tion fjcj) is the identity and f'lf is chain homotopic to the identity in 
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the sense that there is a linear map 9 ■. Vr(G x H) ^ x H) 

satisfying (p'ljj — 'ipcf) = d9 + 9d. 

Corollary 10 (Eilenberg-Zilber). Given two finite simple graphs G, H. 
There is a ehain map ip from x H) = fl*{G) ® Q*{H) to the 

chain complex 12* (G x H) and a chain map (p in the reverse direction 
such that (p and ip are chain homotopic. 

Examples: 

1) Let G be the house graph and H the sun graph S'i^o,o,o- Both have the 
Betti vector (1,1) and the product has the cohomology of the 2-torus 
as it is homotopic to the 2-torus. The graph Gi has dimension 37/24, 
has no(G) = 12 vertices, ni(G) = 18 edges and V 2 {G) = 6 triangles, 
the graph Hi has dimension 1 and vo{H) = 10 vertices and vfiH) = 10 
edges. The product GxH is a graph of dimension 61/24. It has 120 ver¬ 
tices and 480 edges. The kernel of the Laplacian Li on 1-forms of Gi is 
spanned by ( 8 , — 8 , 8 , 8 , —3, —2, —3, 0, —3, —1, —1, —3, —1, 0, — 8 , 2, — 8 ,1), 
the null space of Li on 1 -forms of Hi is spanned by (— 1 , 1 , — 1 , — 1,1 
,1,1,1,0,0). The Laplacian Li on the product graph G x H is a 
480 X 480 matrix. The Laplacian Li on the simplicial set 12^ (G) 0 
12°(77) © 12°(G) © 12^(77) is much smaller as it works on a space whose 
dimension is the set of degree 3 monoids in fg which is vo(G) *vi(H) + 
ni(G) *vo{H). 

2) The cohomology of G 4 has the Betti vector (1,1). The cohomology 
of G 4 X G 4 is (1,2,1). 

3) The cohomology G x is the same as the cohomology of G. Sim¬ 
ilarly, the cohomology of the cylinder T x C obtained by taking the 
products of a couple of circles with the product of a couple of interval 
graphs is the same than of the torus T because the cube, the product 
of interval graphs is contractible. 

4) The Betti vector of G = x ■ ■ ■ x Cr^. with n* > 4 is the k’th row 
of the Pascal triangle (7?(/c, 0), 7?(7,1), B{k, 2),..., B{k, k)). Easier to 
state, Pg{x) = (1 + xfi. 

The reason for the following inequality is that higher dimensional sim- 
plices spawn more vertices in the graph Gi. If the dimension is not 
uniform, then the dimension of the enhanced graph Gi = G x 71i will 
increase. First of all, the construction of Gi allows to dehne a local 
dimension for all simplices x inside G: it is just the usual dimension 
1 + dim(S'(a;)) in the enhanced graph Gi, where the simplices like x are 
vertices. 

Lemma 11 (Dimension inequality). Given a finite simple graph G. 
Then dim(G x Ki) > dim(G). 
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For K 2 , or any two point graph, we have equality. For all three point 
graphs different from the three point graph xy + x + y + z, we have 
equality too, while for xy+x + y + z we have a first inequality dim(G) = 
2/3 and S{G) = dim(G'i) — dim(G) = 1/3. Denote by dim(S'(a;)) the 
dimension of the sphere in Gi and by dim(a;) = (1 + dim(S'(a;))) the 
dimension of the vertex x in Gi. The average of all these numbers is 
dim(Gi). We can not use induction as the unit spheres S{x) are not of 
the form HxKi with a smaller dimensional graph in general. However, 
we can do it with a chain. We therefore prove a stronger statement for 
chains / = ^^e degree 1 monomials be the vertices and the 

sphere S{x) of a vertex x consist of the sum g of degree 2 monomials in 
/ which contain x divided by x. The small dimension dim f f) of a chain 
/ is now recursively dehned as the average of the small dimensions of 
the unit spheres of the vertices. The large dimension dim(/) of the 
chain is the usual dimension of the graph Gf. For a chain / which 
comes from a graph G, the small dimension is the dimension of G, the 
large dimension is the dimension of Gi which is the average over all 
large dimensions of the unit spheres S{x), over all vertices x in Gi. 
The dimension inequality follows from the more general statement: 

Lemma 12 (Generalized dimension inequality). If f is a chain then 
dim(/) > dim ( f). 

Proof. The more general statement can now be proven with induction 
with respect to a lexicographic ordering which looks first at the the 
number of variables appearing in the chain and then at the number of 
monomials. This is a total ordering so that induction works. Unlike for 
graphs, the unit sphere of a vertex x is now a chain with one variable 
less for which we can look at the small and large dimension. Also the 
unit sphere of a simplex a; is a chain for which we can look at the small 
and large dimension. The unit sphere of a simplex has in general the 
same number of variables but less monomials. □ 

Example. 

The chain / = xyzw + xyz + xy + xw + x + y has two vertices x, y. (A) 
Lets look at x first: the unit sphere of x is Sf{x) = yzw + yz + y + w 
which has two vertices y,w. (i) The vertex w in Sf{x) has the unit 
sphere S{w) = yz which has no vertices so that w has dimension 0. 

(ii) The vertex y in Sf{x) has the sphere S{y) = zw z, a chain with 
one vertex which has a 0-dimensional sphere w so that 2 : has dimen¬ 
sion 1, S{y) = z + zw has dimension 1 and y has dimension 2 in /. 
The dimension of the vertex x is therefore (0 -|- 2)/2 = 1. 

(B) The unit sphere of y in / is —xzw — xz which has no vertex so that 
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y has dimension 0. (C) The chain has small dimension (0 + l)/2 = 1/2. 
The two computations (A)-(C) show that dim(/) = 1/2. The large 
dimension of the chain is the average of the large dimensions of the 
monomials in /. It is the dimension of a graph Gf with 6 vertices and 
11 edges. It has dimension dim(Gj) = 2.7333 which is dim(/). 

Remarks: 

1) The inequality looks like the known inequality between Hausdorff 
dimension and inductive dimension in the continuum. But we do not 
see any relation yet. 

2) When iterating the rehnement construction, we get a sequence of 
graphs Gn with a dimension which converges to an integer. 

Lets elaborate on this remark a bit more. The largest simplex will be 
the seed which grows fastest and take over. 

Proposition 13 (Rehnement asymptotics). Start with a finite simple 
graph G = Gq and define recursively Gn = Gn-i x Ki, then dim(G„) —)■ 
d, if Kd+i is the largest cligue in G. 

Proof. The dimensions can only grow or stay by the lemma. Since the 
dimensions are bounded by d for all n, the dimensions must converge 
to some value smaller or equal than d. The limiting value is d because 
the number Vd{Gn) of Kd+i subgraphs of Gn grows exponentially and 
the growth of that large geometric cluster Hn formed by the union of 
the Kd+i simplices takes over: the number vfiGn) of subgraphs 

of Gn is dominated by vfiHn) of iLz+i subgraphs of i7„. □ 

We see that when starting with an arbitrary network, we asymptotically 
get to a geometric graph. This could be interesting if we look at a 
dynamical emergence of graphs through rehnement or evolution. It 
would explain, why geometric and manifold-like structures with integer 
dimension appear naturally: the lower dimensional structures grow 
slower and are overtaken by the growth of largest structures which 
make sure that a larger and larger part of Gn will be geometric. The 
geometric part grows fast like the novo-vacuum in Schild’s ladder where 
physics is described by “Sarumpaet rules” [3] . Egan’s novel starts with 
with the words: In the beginning was a graph, more like diamond than 
graphite. 

The dimension of G x Ki is not necessarily equal to the dimension of 
G. For example, the dimension of the lollipop graph G is 5/2 = 2.5, 
the dimension of Gi = G x iLi is 11/4 = 2.75 is slightly higher. While 
G X Ki has a cover for which the nerve graph is G, the dimensions 
don’t agree because higher dimensional parts of G contribute more to 
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Figure 11. Starting with the house graph G to the 
left, we build Gi = G x K 2 , then G 2 = Gi x Ki etc. 
The graph remains non-geometric but the geometric part 
takes over more and more weight as higher dimensional 
simplices divide into more pieces. In the limit we get to 
a geometric graph. We see that even when starting with 
a random graph, the refinement process moves the graph 
to a geometric graph. This is why geometric graphs are 
so natural. What happens from Gn —t Gn+i is that the 
simplices of Gn become the vertices of Gn+i- As we like 
to think about simplices as basic “points”, this picture 
is a natural evolution. 



Gi. 


We have looked a bit at the statistics and computed the difference 
dim(G X it'i) — dim(G) which is related to the variance of the dimension 
random variable x —)• dimG(a;) on a graph. See Figure (12). 

The preservation of dimension is an essential ingredient of a graph 
homeomorphisms, as it is in the continuum. It is no surprise that in 
the fractal case, we have the same inequality like Hausdorff dimension: 
dim(X X Y) > dim(df) -|- dim(y). Note however that the equality 
dim(G X H) = dim(Gi) -|- dim(iJi) holds in full generality, for any 
hnite simple graph. We use the notation dim(G)(x) = dim(S'G(x) 
where Sg{x) is the unit sphere of the vertex x in G. 


Theorem 14 (Dimension of product). Given two finite simple graphs 
G,H. The dimension of G x H is equal to the sum of the dimensions 
of Gi = GxKi and Hi = H X Ki. More precisely, the dimension of 
any vertex {x, y) in G x H the dimension is the sum of the dimensions 
of X in Gi and y in Hi: 

dim(G X H){x,y) = dim(Gi)(x) -|-dim(Ffi)(|/) . 


Proof. Use induction with respect to the sum of dimension. Fix x G 
Gi,y G Hi and assume dim(S(a;)) = a, dim(il(i/)) = b so that a — 1 = 
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i<5(dim) 

a 



Figure 12. The dimension inequality is illus¬ 
trated quantitatively by computing the discrepancy S = 
dim(G X Ki) — dim(G) for 2000 random Erdos-Renyi 
graphs in E{12, 0.5) in relation to the dimension dim(G). 
The inequality becomes an equality for geometric graphs. 
There, the dimension of G is an integer and agrees with 
the dimension of the enhanced graph Gi = G x Ki. 


dim(S'(x)) and b — 1 = dim(S'(|/)). We only have to show that the 
unit sphere of {x, y) in {G x H) has dimension a -|- 6 — 1 so that (x, y) 
has dimension a -|- 6. The sphere S{x,y) is the union of two graphs 
B{x) X S{y) U S{x) X B{y), which overlap in S{x) x S{y). The two 
graphs have both dimension c = a+b—1 and S{x) xS{y) has dimension 
c — l = a-l-& — 2. Lemma: if two graphs A, B have dimension c and 
their intersection has dimension c — 1 then AU B has dimension c. □ 
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For example, if G is the tadpole graph and H is the house graph, then 
dim(G'i) = 71/44 and dim(G 2 ) = 37/24. The product graph G x H 
has dimension 833/264 = 71/44 + 37/24. 

Corollary 15 (Dimension inequality). IfG,H are finite simple graphs, 
then dim(G x H) > dim(G') + dim(i7). 


Proof. We have dim(Gi) > dim(G) and dim(i7i) > dim(i7). 
result follows from Theorem (14). 


The 

□ 


Given a geometric graph G. A subgraph f/ of G is called convex if 
between any two vertices x, y in U , there exists a shortest connection 
between x, y which is inside U. Note that shortest connections are not 
necessarily unique. The convex hull of a hnite set Y of vertices in G is 
the intersection of all convex subgraphs of G which contain Y. 


Theorem 16 (Topology), a) If G is geometric, then Gi = G x Ki is 
homeomorphic to G in the sense of [2Z] ■ 

b) If G has no triangles, then Gi = G x Ki is homeomorphic to G in 
the classical sense of topological graph theory. 


Proof, a) Assume the graph has dimension k. Let Xi,... ,Xn be the 
vertices of G. We have to hnd a cover {f/(xj)} of Gi such that the nerve 
graph of the cover is the original graph G. As U{xi) we look at the 
unit sphere Scixi) in the original graph and build the smallest convex 
hull in Gi containing all the vertices of this sphere. The intersection of 
two U{x), U{y) is fc-dimensional and nonempty if and only if x,y were 
connected in the original graph. Note that according to the dehnition 
[27] . intersections U{x)P\U{y) which are lower dimensional do not count 
as edges in the nerve graph. 

b) The graph Gi is a rehnement, where every edge has got an additional 
vertex. □ 


Examples. 

1) If G = Gn then the U{xi) are balls of radius 2 in the new graph Gi 
centered at the original vertices. Two sets U{x),U{y) overlap in line 
graphs of diameter 2 if a: and y were connected in G. 


Theorem 17 (Continuity). If G and H are homeomorphic geometric 
graphs and K is a third geometric graph, then G x K is homeomorphic 
to H X K. 

Proof. We use that G and Gi are homeomorphic and that H and Hi 
are homeomorphic. if Ut is a cover of Gi and V) is a cover of Hi such 
that the nerve graphs of are the same and Wj is a cover of Ki, then we 
have a cover (Gj x Wj) and (V) x Wj) of G x K and H x K. □ 
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Remark: We don’t yet know whether this can be proven for arbitrary 
hnite simple graphs. One hrst has to show that in full generality that 
G and Gi are homeomorphic. This is not yet done. 

Let G be a hnite simple graph. Dehne the curvature of a simplex x 
in G as the usual curvature |T6] of the vertex x in the graph G x Ki. 

Theorem 18. The simplex curvature K{x) satisfies the Gauss-Bonnet 
relation '^^K^x) = x(G). The function K{x) is equal to the expec¬ 
tation E[if[x)] when integrating over all colorings f on the simplex 
graph. 

See [IHIEE]. 

Remark: We have already a natural ’’curvature” for a graph located 
on simplices, it is constant 1 on the even dimensional simplices and —1 
on the odd dimensional simplices. But this has little to do with the 
usual curvature. 

Now, we look at symmetries: 

Lemma 19. Given two finite simple graphs G, H. The automorphism 
group Aut(G) is a subgroup of the automorphism group of G x H. 

Proof. Each element T of Aut(G) produces a permutation of the vari¬ 
ables xi,... ,Xn appearing in the algebraic description g of G. If h is the 
algebraic description of H using variables yi,... ,ym, then A acts on fg 
by permutations of the elements xi,... ,Xn and produces a symmetry 
ofGx H. □ 

Remarks: 

1) The case Gi = G x iLi ~ Gg with G = G4 shows that the symmetry 
group can become bigger. 

2) If the vertices of Gi = G x iLi form a group and also the vertices 
of Hi =H X Ki have this property then the vertices of G x H form a 
group, the product group. 

Next we look at homotopy: 

Lemma 20. If G = AUB and H = GUD are two finite simple graphs 
and the subgraph A of G is homotopic to the subgraph G of H and B 
is homotopic to D and Ad B is homotopic to G (1 D, then AU B is 
homotopic to G A D. 

Proposition 21. The graph Gi = G x Ki is homotopic to G. 
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Proof. Use induction with respect to the size of the graph. Its true for 
Ki. Given a graph G with n vertices. Add a new vertex x to get a 
larger graph H with n + 1 vertices. The unit sphere S{x) is in G. Now 
Gi = G X A'l is a subgraph of Hi = H x Ki. Let yi,... ,yk denote 
the simplices in S{x), they are the vertices in Si{x) = S{x) x Ki. The 
graph Hi has new vertices N = {yix,... ,ykx} which were not in Gi. 
Our induction assumption is that Gi is homotopic to G and since S{x) 
is a subgraph of G also Si{x) is homotopic to S{x). But then the 
graph in H generated by S{x) U {t} is homotopic to the graph in Hi 
generated by S'i(a;) U A^. Now H is a union of two graphs S{x) Ux and 
G intersecting in S{x) and Hi is the union of two graphs S'i(a;) U N 
and Gi intersecting in Si{x). Use the lemma. □ 

Theorem 22. Let G, H, K be finite simple graphs. IfG,H are homo¬ 
topic, then G X K and H x K are homotopic. 

Proof. Show it on an algebraic level; if g ^ g + k is a homotopy step, 
then xg —)■ x{g + /c) is a homotopy step on each hbre. □ 

It follows that the product dehnes a group operation on the homotopy 
classes. This monoid dehnes then a Grothendieck group. 

The Gartesian product also dehnes a direct sum on isomorphism classes 
of vector bundles. 

Example. 

1) The product of two trees is contractible. 

2) The product G^ x H is never contractible. 

Finally, lets look at the chromatic number c{G). 

Lemma 23 (Minimal coloring). The graph Gi = G x Kq satisfies 

c{Gi) < c{G) . 

It is minimally colorable. If Kn+i is the largest complete subgraph of 
G, then c{Gi) = n + 1. 

Proof. If a; is a vertex in Gi, denote hj k = Dim(a;) the dimension of the 
corresponding simplex x in G. The function Dim is a coloring because 
two adjacent vertices have diherent dimension. Because Dim takes 
values from 0 to n if K^+i is the largest simplex in G, the chromatic 
number is n + 1. □ 

Remarks. 

1) We see that for graphs Gi we can immediatly compute the chro¬ 
matic number by just looking at the largest clique. 
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2) For geometric graphs G, this means that Gi is minimally colorable 
and therefore is Eulerian in the sense of m For two dimensional 
graphs, this is equivalent to Eulerian in the classical sense. 

We can actually compute the chromatic number of any product graph 
GxH: 

Theorem 24 (Chromatic number of product). IfG has a largest clique 
Kn and H has a largest clique Km, then c{G x H) = n + m — 1. 

Proof. Again just look at the function Dim. It is locally injective and 
is so a coloring. □ 

This result actually has the previous lemma as a corollary because Ki 
has the largest clique Ki so that c{G x Ki) = n + 1 — 1 = n. 

Examples. 

1 The chromatic number of a product of two trees is 3. 

2 The chromatic number of K 2 x K 2 is, ?>. The chromatic number of 
K2 X K2 X K2 is 4. 


6. Examples 

Example. 

Let G be the house graph and H the Lollipop graph. We have dim(G) = 
22/15, 60 (G) = 1, 61 (G) = l, 6 fc(G) = 0,fc > 1 and bo{H) = l,bk{H) = 
0, A; > 0, dim(Lr) = 5/2. The graphs Gi = G x Ki,Hi = H x K 2 
satisfy dim(Gi) = 37/24 and dim(iLi) = 11/4. The graph GxH has 
dimension 103/24. The dimensions of the individual points are given 
by the dimension spectrum dimspec(Gi) ={ 1, 1, 1, 7/4, 2, 2, 2, 7/4, 

1, 2, 1, 2 }, dimspec(i7i) = {3,11/4,3,3,3,3,3,3,1,3,3,3,3,3,3,3,1}. 
The dimension spectrum of G x if is the Cartesian product of these 
two lists: dimspec(G x H) = {15/4, 15/4, 15/4, 4, 4, 4, 4, 4, 4, 4, 4, 

4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 2, 2, 2, 9/2, 19/4, 19/4, 19/4, 19/4, 

5, 5, 5, 19/4, 5, 5, 5, 19/4, 5, 5, 5, 19/4, 5, 5, 5, 19/4, 5, 5, 5, 19/4, 

5, 5, 5, 19/4, 5, 5, 5, 11/4, 3, 3, 3, 9/2, 15/4, 19/4, 19/4, 4, 5, 19/4, 

4, 5, 19/4, 4, 5, 19/4, 4, 5, 19/4, 4, 5, 19/4, 4, 5, 19/4, 4, 5, 11/4, 2, 

3, 15/4, 4, 4, 4, 4, 4, 4, 4, 2, 19/4, 5, 5, 5, 5, 5, 5, 5, 3, 4, 4, 4, 4, 4, 

4, 4, 4, 4, 4, 4, 4, 19/4, 5, 5, 5, 19/4, 5, 5, 5, 19/4, 5, 5, 5, 19/4, 5, 5, 

5, 19/4, 4, 5, 19/4, 4, 5, 19/4, 4, 5, 19/4, 4, 5, 4, 4, 4, 4, 5, 5, 5, 5, 4, 

4, 4, 4, 4, 4, 19/4, 5, 5, 5, 19/4, 5, 5, 5, 19/4, 4, 5, 19/4, 4, 5, 4, 4, 5, 

5, 4, 4, 4, 19/4, 5, 5, 5, 19/4, 4, 5, 4, 5, 2, 2, 2, 11/4, 3, 3, 3, 11/4, 

2, 3, 2, 3}. This list had been computed directly by computing the 
dimensions of each unit sphere in G x ii. The Kiinneth formula tells 
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Figure 13. x Ki, x Ki) x Ki, x Ki) x 

K,) X K,. 


Example. 

1) For an octahedron G, the ring element is /g is a + 6 + a6 + c + ac + 
be + abc + d + ad + bd + abd + e + ae + ce +ace + de + ade + / + &/ + c/ + 
bef + df + bdf + ef + cef + def. The graph G x iFi is a barycentric 
subdivision of G of the same dimension. For a geometric graph G, the 
graph G X Ki is geometric again. 

2) Let O be the octahedron, the product O x G 4 is a 3-dimensional 

graph. It is a discrete incarnation of the 3-manifold x S^. 
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Figure 14. The graph Lq x Lg where Ln is the line 
graph. This is a discrete square. The second picture 
shows L 3 X L 3 X L 3 , a discrete cube. 



Figure 15. The lollipop graph G is partly 2- 
diniensional, partly 1 -dimensional so that it is a graph 
with fractal dimension. The graph G x L 3 is partly 2- 
dimensional and partly 3-dimensional, also a fractal di¬ 
mension. 

3) A triangle is given by it's = x + y + z + xy + yz + zx + xyz. Lets 
take the product with K 2 represented hy u + v + uv. It becomes a 
3-dimensional graph. 

4) The product of with K 2 is a 3-dimensional graph. It is a 3- 
dimensional ball of radius 1 with one interior point and all other points 
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Figure 16. The product O x Ki where O is the octa¬ 
hedron graph. The second picture shows I x Ki, where 
I is the icosahedron graph. 



Figure 17. The product O x where O is the octa¬ 
hedron graph. This is a discrete 3 manifold. The second 
hgure shows the product of the sun graph S' 2 , 2 , 2,2 (a cyclic 
graph (74 with rays) and with K 2 . 

on the boundary. The curvature of the interior point is 0, the curvature 
of nine of the boundary points are 1/6 and three are —1/6. By Gauss- 
Bonnet, the sum of the curvatures is 1, which is the Euler characteristic 
of the ball. 



48 


OLIVER KNILL 




Figure 18. To the left, the product of a wheel graph 
with a line graph: both factors are homotopic to a point 
so that the product is homotopic to a point. To the 
right, the product of a sun graph S'(4,4,4,4) with a cyclic 
graph C^. both factors are homotopic to a circle so that 
the product is homotopic to a torus. 




Figure 19. The product of two star graphs Fs is a 
2-dimensional graph. The product of (^4 with the star 
graph 5*3 is 2-dimensional. It is not geometric even so it 
looks like a cobordism between C 4 U C 4 and C^. 
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Figure 20. The product of with K 2 is a 3- 
dimensional ball of radius 1. Its boundary is a 2- 
dimensional sphere. The second picture shows the prod¬ 
uct K 3 with K 3 . It is a 4-dimensional contractible graph 
whose boundary is a 3-dimensional sphere. 




Figure 21. The product of with C 4 is a discrete 
2-dimensional torus. The traditional cartesian product 
C 12 with C 5 is a 1 -dimensional graph as it has no trian¬ 
gles. 





50 


OLIVER KNILL 




Figure 22. The graph product of the house graph G 
with H = K 2 . We see hrst the enhanced graph Gi = 
G X Ki and then the graph G x H. 
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7. Lose ends 

Dimension: 

What is the expectation of X{G) = dim(Gi) — dim(G) on the Erdos- 
Renyi probability space. This random variable measures in some sense 
the distance from a geometric graph as for geometric graphs, the value 
is zero. The value of X [G) is related to the variance of the dimension 
spectrum, dim(x) which is a random variable on vertices of G. It would 
be nice to quantify this more. 

Hausdorff dimension: 

Due to the analogy of results with Hausdorff dimension dimj^, one can 
ask whether there is more to it and whether for any subsets X, Y of 
Euclidean space R^, there are sequences of graphs Gn,Hn such that 
dim(G„) —)■ dim//(W) and dim(i7„) —)■ dim/^CE). 

Curvature: 

The construction G ^ G x Ki gives a convenient natural curvature 
on the simplices of a graph G. Gauss-Bonnet holds as they add up to 
Euler characteristic. It is also true that curvature is the average over 
all functions [181 ES] 

What is the relation between the curvature of G, H and G x HI The 
example of x Gm shows that the product of two zero curvature 
spaces can develop curvature. The case (W 3 x Ki) x Ki shows that the 
curvature can become very negative. 

Spectrum: 

What is the relation between the spectrum oiGxH with the spectrum 
of G and HI We know that the Laplacians live on different spaces. The 
eigenvalues appear to correlated to the eigenvalues A/i, where A is an 
eigenvalue oi L{G x Ki) and /i is an eigenvalue of L{H x Ki). 

Chromatic number: 

The case G^ —)■ Cio shows that the chromatic number can become 
smaller. The product of two geometric graphs is Eulerian in the sense 
[29| which shows that if G is a geometric graph, then Gi = G x Ki 
is Eulerian. The map G —)■ Gi regularizes in some way as for geo¬ 
metric graphs, the unit spheres of Gi are Eulerian spheres. Of course, 
it would be nice to see a relation between the chromatic polynomials 
c(G) and c(Gi) of G and Gi. It is not always true that cg{x) divides 
cgi(t) but it is often the case. For the graph G = for example, 
cg{x) = — 3x‘^ + 2x and cgj (t) = 62*x — 191 * x^ + 240 *x^ — 160 * 
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+ 60 * X® — 12 * X® + x'^. Now cgi(x)/cg(x) = 31 — 49 * x + 31 * 
x^ — 9 * x^ + x^. For G = K4, this chromatic fraction of the graph G 
is -8834963 + 21070773x - 23177760x2 + 15536807x3 - 7055131x^ + 
2278645x3 - 533896x3 ^ 90682x^ - 10932x8 + 890x® - + x^h Al¬ 

ready for G = G 4 , the chromatic fraction is not a polynomial. Here 
is more random example, where the chromatic fraction is a rational 
function and not a polynomial is the graph G with 7 vertices and edges 
{(1,4), (1, 6), (2,4), (2, 6), (2, 7), (3,6), (4, 6), (5, 7), (5, 8), (6, 7))} for which 
the chromatic fraction is only rational. Still, the chromatic number is 
3 for both graphs. With respect to the known relation c{Gi) < c{G), 
we have inequality for all graphs G 2 n+i- For complete graphs Kn, both 
sides are n-|-1, for trees, both sides 2. For geometric graphs, all graphs 
Gi have the property that c{Gi x Kq) = c{Gi). On the list of all 38 
connected graphs with 4 vertices, we have equality c(Gi) = c(G), on 
the list of all 728 connected graphs with 5 vertices there are a dozen for 
which the chromatic number has dropped, the reason always being a 5- 
cycle of chromatic number 3 becoming 2-colorable after the rehnement. 

Constructing spheres: 

What kind of spheres appear as intersections of unit spheres in discrete 
tori Gni X • • • X Gni, with nk > 4. We have seen that they all are unions 
of solid cylinders glued along a torus. Assume we have obtained a class 
of spheres like that. Now take the product of such spheres and take 
again unit spheres. We get a larger and larger class of spheres. Can 
one characterize them? Similarly, we can ask what kind of spheres can 
be obtained by applying the join operation to spheres, starting with a 
smaller class of spheres. 

Determinant: 

Can we say something about the complexity Det(L) of the product, 
where Det is the pseudo determinant [31]. If G = = Ki, then 

the complexities of Gi, Hi are 716800 and 3. The complexity of Gi x G 2 
is already 99446661202899312347583662774063513018614477321403222712572446349721600. 
The pseudo determinant of the Laplacian of Gi is 2240 and of Hi is 3. 
The pseudo determinant of the Laplacian of Gx if is 21101889958483152. 

Homotopy classes: 

If Gi is homotopic to G 2 and Hi is homotopic to H 2 , then Gi x G 2 
is homotopic to Hi x H 2 . The sum and product therefore produces a 
sum and product in the ring of homotopy classes. The contractible ring 
elements form a class representing the 1 element. The empty graph is 
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Figure 23. A graph G and the enhanced graph 
Gi = G X Kq. Both graphs have chromatic nnmber 3. 
The chromatic polynomial of G is (—2 + 1 + x)'^x, 

the chromatic polynomial of Gi is (—2+a;)(—l+a;)^a;(7 — 
5x + a;2) (12193 - 43730a; + 72229x^ - 72238a;3 + 48366a;^ - 
22614a;® + 7465a;® - 1715a;^ + 262x® - 24a:® + a:^®). The 
chromatic fraction is not a polynomial. Still, the chro¬ 
matic nnmber is 3 for both graphs. 


the 0 element. Since also SGi is homotopic to 6 G 2 - Homotopy ex¬ 
tends to the ring of chains. Dehne contractibility in the same way by 
telling G is contractible if there exists x snch that both S{x) and G{x) 
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are contractible where G{x) is the non-constant part of the polyno¬ 
mial where x is replaced by 1. The unit sphere is the unit sphere of 
the graph Gf. Homotopy is so dehned algebraically as reduction steps 
f ^ f — S{x) with contractible S{x) or reverse steps. Contractibility 
is dehned for / in a polynomial ring as the property that there is a 
variable x* such that S{xi) and / —S'(xj) are both contractible. For ex¬ 
ample a line graph xy + yz + x + y + z is contractible as S'(z) = z + yz 
and f — S{z) = xy + x + y are both contractible. The cycle graph 
xy + yz + zw + wx + x + y + z + w is not contractible. Also all its unit 
spheres like S{x) = xy + wx or S{y) = xy + yz which represent two 
point graphs are not contractible. 



Figure 24. The product of two trees is a 2-dimensional 
graph which is contractible as it is the product of two 
graphs which are homotopic to a point. 

Suspension and join; 

The topological notions of “suspension” and “join” needs the Cartesian 
product. The dehnition of a suspension of a graph G is SG = G x 
//{(x,0) ~ (|/,0 ),(x, 1) {y, 1 )}. Also in the graph case, it produces 
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from a /c-sphere a (/c + l)-sphere. More generally, the join of two graphs 
G, H can be defined as 

{Gx H X K2)/{{x, yi, 0) ~ {x, y2, 0), (xi, y, 1) ~ (xs, yG)} ■ 

An obvious questions is whether results like the double suspension the¬ 
orem of Cannon-Edwards can be proven directly in graph theory. It 
certainly is true, as one can see by building corresponding manifolds 
from the graphs. But it would be nice to have a direct discrete proof. 
The double suspension of a graph G is S'^G = (GaS'o)aS'o. A homology 
n-sphere is a graph G which is geometric of dimension n, which has 
the Poincare polynomial Pg{x) = 1 + x^ but which is not a n-sphere. 
The discrete Cannon-Edwards theorem tells that if G is a homology 
sphere, then H = S'^G is a sphere. To prove this within graph the¬ 
ory, one has to show that every unit sphere in if is a (n — l)-sphere 
and that removing one vertex in H produces a n-ball, a contractible 
n-dimensional graph with n — 1 dimensional sphere boundary. 

Higher order curvature; 

It is natural to inquire whether for 2-dimensional geometric graphs 
Gi = G X ifi, the higher order curvature K{x) = 2Si{x) — 5*2 (t) pro¬ 
duces a Gauss-Bonnet result. This is not the case as for the octahedron 
G, the graph Gi has a sum of the second order curvatures which is 0. 
The culprit are the degree 4 vertices. 

The sequence of graphs G„. 

What happens with G„ dehned recursivly by G^ = G„_i x Ki is 
that the minimal curvature goes to — cx), at least if Go is two di¬ 
mensional. It would be interesting to know how the curvature spec¬ 
trum ^(Gn) = [min(A'G„(a^), Hiax(A'G„(a^)] grows for higher dimensional 
graphs. The numerics is tough as the graphs sizes explode exponen¬ 
tially. L If Go = ffs, then the curvature spectra are k(Go) = [1/3,1/3], 
a(Gi) = [0,1/6], a(G 2 ) = [—1,1/3],.... If Go = ff 4 , then k{Gq) = 
[l/4,l/4],«:(Gi) = [0,1/6]. 

Global De Rham: 

We have constructed a chain homotopy between de Rham and simpli- 
cial chains. If we dehne a discrete d-manifold as a graph which locally 
can be wrtten as a product Ui = Gn x ■ ■ ■ x Gid of networks such that 
on the intersection Ui fl Uj there are coordinate change graph home- 
omorphisms, then the chain homotopy can be pushed through. We 
don’t even need the graphs Gij to have to be 1-dimensional, nor have 
they be to be geometric. The combinatorial de Rham theorem allows 



56 


OLIVER KNILL 



Figure 25. In this figure, we have taken joins by tak¬ 
ing the interval K 2 . The join Cg * S'o of a cycle graph 
with a O-dimensional sphere is a 2-diniensional sphere 
seen to the left. The join O* 5*4 of the octahedron with a 
circular graph is a 4-diniensional geometric sphere seen 
to the middle. The join O * O of the octahedron graph 
O with itself is a 5-dimensional geometric sphere. It has 
no = 728 vertices, vi = 14168 edges, 72960 triangles. 


to deal with such discrete manifolds more effectively similarly as the 
continuum de Rham theorem does. 

Moving frames: 

While a full classical Hopf-Rynov theorem on graphs is impossible be¬ 
cause of the finiteness of any reasonable tangent space and therefore, 
a quantum flow is needed for a reasonable notion of geodesic map. 
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Figure 26. This figure shows more joins. This time, 
we replace K 2 with a longer line graph representing 
an interval. The first case is Cs a Sq with interval graph 
L4. It is a 2-dimensional sphere. The second case is 
K 2 ^K 2 with interval graph L4. It is a three dimensional 
ball with two dimensional sphere as boundary. The third 
example is a product of two Erdos-Renyi graphs with 5 
vertices, edge probability 0.2 and using the interval graph 
L3. 

there is weaker discrete Hopf-Rynov theorem which tells that under 
some conditions, there is a global unique flow which locally minimizes 
length [28] . What is needed for a geometric graph of dimension d is an 
Eulerian condition which is equivalent to the graph being minimally 
colorable with d + 1 colors. In |28|, we made a stronger assumption 
and required all unit spheres to be projective in order to associate to 
an incoming ray an outgoing ray. There is an other approach which 
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implements a “moving frame” idea of Cartan and works for all Euler- 
ian graphs and in particular with Gi = G x Kq of G x if if G, hf are 
geometric. It in particular does not need the projective assumption. 
A discrete moving frame is a pair where a is a d-dimensional 

simplex containing x. The frame propagation dehnes a geodesic flow 
in an Eulerian graph without the need of a projective involution: it 
is described inductively with respect to dimension: use the geodesic 
flow on the Eulerian unit sphere S{x) to flip the frame {x^a) on the 
sphere S{x) to the anitpodal side (x'jcr'), where either the simplex a' 
or the vertex x' has maximal distance (both cases are possible as the 
example in Figure (27) shows), then reflect the vertex x' along the face 
a' \ {x'} to a vertex x" on the adjacent simplex a". This propagation 
(x, a) —)■ (x", a") dehnes a variant of a geodesic how for which the 
simplex a plays the role of the direction, as well (if we keep track of 
the frame in the simplex) produces a parallel transport on the graph. 
The point we want to make is that for all product graphs G x if and 
more generally for all discrete manifolds obtained by patching product 
graphs, where each factor is a geometric graph, there is a canonical 
notion of geodesic how and parallel transport. In some sense, there is 
a unique Levi-Civita connection on such graphs. 


Figure 27. The stellated cube G is a sphere with 14 
vertices. It is Eulerian. The graph Gi is also Euler¬ 
ian with 74 vertices. While G has a natural antipodal 
map (there is a unique opposite vertex with maximal 
distance), this is no more the case for Gi. The antipo¬ 
dal point of some vertices is a triangle. The discrete 
analogue of a moving fame associates to a vertex in the 
d-dimensional graph a d dimensional simplex. 
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